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b) Resistive heating of the field-emission centres on the cathode 
leads to a rupture of one or more of these projections 3~’). 

c) Positive ions, produced at the anode by field-emission elec- 
trons, locally heat the cathode ”)8). 

d) Loosely attached pieces of material (clumps) on either elec- 
trode, carrying a charge, are torn from the electrode surface and 
cross the gap under influence of the field. On striking the other 
electrode they may cause high temperatures and local evap- 
oration 9). 

e) Breakdown occurs, when at a certain voltage a free electron 
produces an avalanche of charged particles. Primary electrons may, 


(3) 


a] 
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Fig. 1. Breakdown voltage V, as a function of gap length d. 

1. Trump, J. G. and R. J. Van de Graaff, 12), iron electrodes, impulse 
voltage. 

2. Mason, R. C.}5), tungsten cathode, copper anode, direct voltage. 

3. Leader, D. 18), iron electrodes ,direct voltage. 

4.. Rosanova, N. B. and V..L. Granovskii 17), nickel cathode, molyb- 
denum anode, impulse voltage. 

5. Rosanova, N. B. and V. L. Granovskii 17), molybdenum electrodes 
alternating voltage. ; 

6. Denholm, A. S. 7), iron electrodes, alternating voltage. 

7. Wijker, W. J., a) iron electrodes, b) copper electrodes, alternating 
voltage. 

8. se W. J., stainless steel electrodes, impulse voltage (d = 0.1 mm 
only). 
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produce positive ions and eventually photons at the anode; these, 
in their turn, striking the cathode, may produce electrons and so 
the process continues 10-18), 

Because of the great number of variables in the problem, no 
theory up to now describes quantitatively all experimental results. 

Some measurements of the breakdown voltage as a function of 
gap width are represented in fig. 1. The slope of the curves varies 
between 0.5 (as predicted by Cranberg’s theory 9)) and 1. The 
voltage gradient in the gap decreases from about 2 MV/cm at a gap 
width of 0.1 mm to about 60 kV/cm at 10cm. The spread in the 
measurements may be caused by differences in degassing methods, 
in the material and in the shape of the electrodes, as well as in the 
electrical circuits used (direct, alternating and impulse voltage 
circuits). 

In most previous investigations less attention has been given to 
the pre-breakdown phenomena. It is to be expected, however, that 
the pre-breakdown current in the gap is important for the final 
breakdown. Measurements of this current will be described in this 
article, both for impulse voltage and for 50 Hz alternating voltage 
breakdown. §2 deals with the experimental arrangement. The 
conditioning of the electrodes will be described in § 3. The results of 
the measurements are dealt with in the §§ 4 and 5, while in $6 a 
discussion of the results is given. 


§ 2. Experimental arrangement. The experiments have been 
carried out in two different vacuum chambers, one being a de- 
mountable, rubber-sealed vessel, evacuated by a mercury diffusion 
pump, the other being an all-glass system, directly sealed to an oil 
diffusion pump. In both cases a liquid-air trap was placed between 
the pump and the vacuum chamber. With the available pumps a 
working pressure of 6 x 10-5 and 8 x 10-6 mm Hg respectively 
could be established. This corresponds to a mean free path of at 
least 80 cm for molecules (air) and 500 cm for electrons, which is 
large compared with the dimensions of the gap. 

In the demountable apparatus, the distance between the elec- 
trodes could be varied during the experiments by means of a rubber 
membrane, while in the all-glass system this was done by means of 
an external magnet. Figs. 2 and 3 give cross-sectional views of the 
vessels. The electrodes in each vessel were of a similar shape and 
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replaceable. The diameter of the electrodes in the demountable 
apparatus was 26mm and the convex surfaces had a radius of 
curvature of 40 mm. For the all-glass system these figures were 18 
and 25 mm respectively. Copper, iron and stainless steel electrodes 
were used. 


PUMPS 


ee ee 
0 50mm 


Fig. 2. Cross-sectional view of the demountable vacuum chamber. 


The applied tension was either a 50 Hz alternating voltage, 
adjustable from 0 to 45 kV peak value, or an impulse voltage with 
a peak value of about 32 kV and a rise time varying from 1.0 us 
tot 1.3 ms. In the a.c. circuit a current limiting resistor was placed 
in series with the gap. The gap voltage was measured with a 
resistive voltage divider in parallel with the gap. The current was 
recorded with the aid of a 1000 Q resistive element. 

The impulse voltage was generated by charging a 120 pF con- 
densor through a rise-time defining resistor of variable value 
(Rz, fig. 4). The charging of this condensor was started by triggering 
the spark gap SGI, which is initiated by the discharge of a con- 
densor (2 uF) through a thyratron (EC 50) and the primary of an 
induction coil. This caused the small spark gap SG II to fire. The 
middle sphere of SG I then fell to earth potential and SG I broke down. 
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VACUUM GAP 


SUPPLY 0.04 pF CRO. 


| 235M 
—T/ME BASE 
[] 10m 


TIME BASE 
Fig. 4. Impulse generator with starting circuit (see text). 


6 W. J. WIJKER 


The thyratron was fired at a time interval after a starting im- 
pulse had been given to the time base of the cathode ray oscilloscope 
(c.r.o.) by switching S from load to start position. This time interval 
could be adjusted by the condensor Cy. However, it had a minimum 
value, owing to the self-inductance of the induction, coil. Therefore, 
for small rise times (< 10 us), the time base was started by a pulse 
taken from a resistor (10 MQ) in series with the voltage divider on 
the supply condensor (0.04 uF). The delay then could be adjusted 
by the small variable high-voltage condensor Cy2. 

The gap voltage was measured with a capacitive voltage divider, 
the largest condensor of which was specially constructed to have a 
negligible self-inductance in order to avoid false readings at small 
risetimes. 

Current was measured by means of a low-inductance 10 kQ 
resistor. All measuring voltages were photographically recorded 
with a c.r.o. 


§ 3. Conditioning of the electrodes. Before mounting, the elec- 
trodes were polished, washed in benzene and then in destilled water. 
Care was taken not to touch the electrode surfaces after cleaning. 

Outgassing of the electrodes was done by induction heating to a 
temperature of well over 700°C for 5 to 6 hours. In the all-glass 
system the electrodes could be heated in situ, but the electrodes of 
the demountable vessel had to be outgassed outside the vessel. This 
was done in special glass containers connected to the pumping 
system. After heat treatment, these containers were broken and 
the electrodes mounted in the vacuum chamber. The experiments 
showed that exposing outgassed parts to the atmosphere for a 
short time (not exceeding a few hours) had no influence on the 
results of the measurements (see also 14)). If a voltage was applied 
to the electrodes heated outside the discharge chamber, a kind of 
gaseous discharge would first appear. After lowering the voltage 
and increasing it again, only spark-like discharges could be observed 
at higher values of the voltage. Provided the pressure did not rise 
above 0.01 mm for long periods, the gaseous discharges did not re- 
appear. With the electrodes that did not come in contact with the 
atmosphere after heating no such discharges were observed. 

Degassing of the electrodes by means of a low-current hydrogen 
discharge (5 mA) at pressures from 0.1 to 1.0 mm Hg has also been 
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tried (gap setting 3mm). Preliminary a.c. experiments showed, 
however, that increasing the voltage across the gap always resulted 
in a gaseous discharge of high intensity, either in the space between 
the electrodes or in the whole vessel. These discharges caused 
the pressure in the vessel to increase considerably, sometimes to 
as much as a hundredfold of the initial pressure. Better results 
might be expected if the current in the hydrogen discharge were 
much larger (order of 1 A; see *)). 

For the first few voltage runs, small sparks were always observed 
between freshly mounted well-degassed electrodes, at relatively 
small values of the voltage (about 10 kV at 2 mm gap width). These 
sparks had the appearance of tiny red to white hot particles crossing 
the gap and hardly affected the gap voltage (a.c. experiments). 
These are probably small pieces of electrode material, weakly 
attached to the surface as a result of the previous machining, that 
are torn loose by electrostatic forces 9). Later on, these sparks no 
longer appeared with the same set of electrodes, except sometimes 
in the voltage interval just prior to breakdown (see § 5). 


§ 4. Impulse voltage breakdown (all-glass system, stainless steel 
electrodes). 

4.1. Breakdown voltage. The effect of the following items 
on the breakdown strength of a 0.1 mm gap has been investigated: 
a) preceding breakdowns, 

b) the time interval between the breakdown under observation 
and the foregoing one, 
c) the rise time of the impulse voltage. 

4.1.a. A graph of the breakdown voltage Vy versus the number 
of breakdowns after the mounting of a particular set of electrodes 
is shown in fig. 5. This graph holds for a series of breakdowns with 
a time interval between the discharges of 30s. This interval also 
has an effect on the breakdown voltage (see 4.1.0). Therefore the 
number of breakdowns, required to reach the maximum value of 
the breakdown voltage again, depends on the time elapsed since 
the last breakdown occurred and on the total number of break- 
downs made on the same set of electrodes. Thus, for much used 
electrodes that had not been experimented with for a few days the 
maximum level was reached after about 10 breakdowns. For 
relatively new electrodes this number is much larger (fig. 5). 
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It may be noted that, once a few tens of discharges had been 
made, the low initial value of the breakdown strength of new 
electrodes never occurred again with the same set. Apparently the 
conditioning of the electrodes by discharges has a permanent and a 
temporary component, the temporary conditioning of much used 
electrodes being achieved with fewer discharges than that of less 
used electrodes. 


Vpin kV 
30 


AVERAGE VALUE FOR WELL COND 


IONED ELECTRODES , 


25 


20 


0 100 200 300 400 500 
NUMBER OF PREVIOUS BREAKDOWNS 


Fig. 5. Breakdown voltage V, as a function of the number of previous 
discharges. Iron electrodes, d = 0.1 mm. Rise time to breakdown 5 lOn ase 


4.1.6. Some experiments have been performed to estimate the 
time in which the temporary conditioning is lost while maintaining 
vacuum conditions. A gap between two well-conditioned electrodes 
showed a decrease in the breakdown voltage if the time interval 
between two successive discharges was of the order of 1 minute or 
longer and the rise time of the voltage was less than, say, 40 us 
(table I). With a larger rise time, the time interval between suc- 
cessive breakdowns could be increased to between 5 and 15 minutes 
before a decrease in the breakdown voltage was obtained (table II). 
Time intervals of more than 5 hours caused a further decrease to 
about the minimum breakdown voltage found with much used 
electrodes. , 

Table I shows a comparison between maximum breakdown 
voltage and breakdown voltage after 5 minutes for different values 
of the rise time (average of 3 observations for each rise time) 
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TABLE I 
A B c 
(us) (kV) (kV) 
2 23 + 0.8 20 + 0.9 Some influence of 5 
5 19 + 2.2 16 + 1.0 min time interval 
20 24 + 1.7 21 + 0.6 
40 26 + 1.0 26 + 0.9 No influence of 5 
85 27 + 0.8 27.5 + 0.5 min time interval 
220 27 +.0.8 26.5 + 0.7 | 


A = rise time to breakdown. 

B = maximum breakdown voltage (11th discharge of a series with an interval of 15s 
between successive discharges). 

C = breakdown voltage 5 minutes after B. 


TABLE II 
A B C 

i 5 min | 22.5kV 26.7 + 1.0kV 
2 15,, SFG ier 25.24 0.4.5, 
3 30,, ick = 25> 4. 0.9%, 
4 60,, 24.0 ,, 26 NOT. 3 
5 2h oo ea 25 Ra 2 
6 Bt 20.5 ,, 26.8 4/05 4; 
7 A 21S... 25.4 det 
8 | oe 22:0 7%, 25. AOs5 
9 ae geo 25 hii 
10 7. Weems 25 3 10.9 


A = time interval after the preceding series (having maximum values of the breakdown 
voltage). 

B = breakdown strength of the first discharge in a series. 

( mean breakdown strength of the following 5 discharges (caused within 5 minutes). 


I| 


Table II, column B, shows the reduction that occurs in the break- 
down voltage of the first breakdown after the time intervals men- 
tioned in column A (rise time to breakdown 1.2 ms). From column C 
it appears that the average value of the breakdown voltage in the 
next five discharges is again equal to the maximum attainable 
value. This indicates a rapid re-conditioning of the electrodes after 
a temporary deterioration. 

It has to be added that if the discharges are forced to follow each 
other with time intervals of less than about 1 ms, then the first of 
such a series generally shows the highest breakdown voltage. This 
is shown by six series of rapidly following discharges with the rise 

‘time of the voltage up to 20kV equal to 0.48 ms. The average 
breakdown voltage of the first discharges of each series was 22.3 + 
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0.6kV, whilst the average breakdown strength of 32 following 
discharges (taken from all six series) was 18.5 + 3.1 kV. For 6 
series with the rise time of the voltage to 20 kV equal to 0.13 ms 
these figures were 24.3 + 0.8kV and 20.9 + 1.7 kV (average of 41 
discharges) respectively. Such a rapid succession can be forced by 
placing a suitable resistor in series with the gap. 


Vp in kV 


0) = | 


X 50 cps 


=2 


iO, fa= jou 10> 10-75 


Fig. 6. Breakdown voltage Vy versus the rise time of the voltage. The 
breakdown voltage for 50 Hz. alternating voltage is indicated by (X). 
Stainless steel electrodes. d = 0.1 mm. 


4.1.c. The relation between breakdown voltage and rise time is 
shown in fig. 6, where three series of observations with different 
rise times are plotted. To make sure that the conditioning of the 
electrodes would be approximately the same for the different series 
(taken on different days), 50 discharges were caused just before the 
beginning of each series. The breakdown voltages in the series were 
each taken from the last one of a sequence of ten breakdowns, 
closely following each other (invervals of 30 s). 

The average value of a few observations with slowly increasing 
direct voltage (1 kV/s) and the maximum breakdown voltage with 
50 Hz alternating voltage are indicated on the same graph. The 
breakdown voltage initially decreases with increasing rise time to 
a.minimum at about 5 x 10-®s; it then rises until a maximum 
occurs between 5 x 10-5 and 1 x 10-4 s; after that it slowly falls 


~~ rr 
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again. For the shortest rise times and those longer than 5 x 10-5s, 


the spread in the measuring series is considerably less than that 
around the minimum. 


4.2. Pre-breakdown current. The pre-breakdown current 
consists of a leading capacitive component and a component in 
phase with the applied voltage. The latter can be found by sub- 
tracting the (calculated) capacitive component from the total 
current. This apparently leads to inaccuracies at the beginning of 
the appearance of the pre-breakdown currents, where the capacitive 
current has its highest value. 

As table III shows, the current intensity, just before breakdown 
is initiated, increases with the number of breakdowns as long as the 
maximum level in the breakdown voltage has not yet been reached. 
Once the maximum breakdown voltage has been reached, this 
current remains of the same order; variations of 20 to 30% may 
occur. 


TABLE. Lil 


(Rise time: O—30 kV in 250 us) 


A | B 6 
1 135 pA 19.5 kV 
2 150.5, EON ok 
3 150 ,, 23.26 
4 : 180 ,, a On 
5 200 ,, | 220 ee 
10 600 ,, 26.8 ,, 
19 | 240 ,, ZO 
20 450 ,, 25.9 ,, 
21 600 ,, 26:8, 
36 1200 ,, i fee 
37 1050", Oey 
38 1100 ,, | 272. 


A = number of discharges 
B = current just prior to breakdown 
C = breakdown voltage 


For higher numbers of discharges the current remains in the 
order of 1 mA. The current-voltage relation prior to breakdown for 
different rise times is given in fig. 7. Because of the large capacitive 
component of the observed current, the values of the curves for the 


two smallest rise times are not very reliable, except for those parts 
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Fig. 7. Current-voltage relations for different rise times of the impulse 
voltage (indicated in brackets). Stainless steel electrodes, d = 0.1 mm. 


10-* 3 107° 3 107* 


eee 10? 
risetime in s 


Fig. 8. Current just before breakdown as a function of the rise time of the 
voltage. Stainless steel electrodes, d = 0.1 mm. 
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just prior to breakdown. Fora rise time of 1 us the current increases 
to 30 mA before breakdown occurs. With increasing rise time, the 
current intensity just before breakdown decreases inversely 
proportional to y/# (fig. 8). The energy flux into the electrodes, 
calculated from the observations given in fig. 8, is in first approxi- 
mation proportional to 1/? (assuming that the fraction of the total 
energy which penetrates into the electrodes is independent of the 
voltage). The heat flow from a source of high temperature, situated 
on the electrode surtace, into its volume is also proportional to 
/t (ref. 18), p. 91). Hence, it seems that a critical energy flux to 
the anode is needed to initiate the discharge. 

Current observations showed that sometimes the current would 
increase when in the tail of the voltage wave the voltage was 
already decreasing. While in some cases breakdown followed, this 
was not always the case. The final transition from pre-breakdown 
current to discharge current (not visible on most photographs) 
happened in a very short time that is estimated to be less than 
LOS. 

The appearance of the electrode surfaces after many breakdowns 
showed that a large amount of anode material was deposited on the 
cathode (fig. 9). Only a negligible amount of material left the gap by 
radial diffusion. A minor part of the anode material appeared as 
droplets on the cathode surface which were rather weakly attached 
and easily removed. From the colours of the electrodes it appeared 
that the whole of the anode had reached far higher temperatures 
than the cathode. 


4.3. Conclusions. The experiments described in §§ 4.1 and 4.2 
have shown that the electrode surfaces are best conditioned by a 
series of breakdowns. Apparently these surfaces are changed by the 
evaporation of projections. This lowers the effective field strength 
at the electrode surface and decreases the effective absorbing 
surface, thereby affecting any surface layers. This leads to permanent 
conditioning. Temporary conditioning is caused by removing the 
surface layers. 

These are re-formed after a time which can be deduced from 
table II. Obviously the removal of the surface layers also takes place 
before breakdown. This makes it understandable that for the 
shorter rise times less time is needed to lower the breakdown 
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voltage than for the larger rise times. The energy needed to remove 
the surface layers, delivered before breakdown takes place, will be 
small at short rise times, thus only allowing the removal of an 
incomplete (formed in a short time) surface layer. 

Analogous reasons can explain the dependence of the breakdown 
voltage on the rise time of the voltage. For very short rise times the 
breakdown voltage can be expected to be high because, again, only 
a small amount of energy is produced at the electrode surface. For a 


Fig. 9. Appearance of the electrodes after a large number of discharges 
(impulse voltage). Left: anode. 


rise time of 1 us the pre-breakdown current must grow to about 
30 mA in order that sufficient anode material is evaporated for 
discharge initiation. Thus, for increasing rise times one would ex- 
pect the breakdown voltage to fall. It is thought however, that a 
second process runs parallel with the first mentioned, this second 
process being the pre-breakdown removal of surface layers that tend 
to decrease the breakdown voltage. The pre-breakdown removal of 
surface layers is thought to commence in the region of 5 us. This 
initially results in an increase of the breakdown voltage with longer 
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rise times. If the rise time is larger than about 100 us, the cleansing 
effect is no longer affected by the rise time and the breakdown 
voltage gradually decreases to its final (d.c.) value. 


§ 5. Measurements with 50 Hz alternating voltage. 

5.1. Breakdown voltage. The breakdown voltage as a func- 
tion of gap width for well-conditioned iron and copper electrodes is 
given in fig. 1, curve 7a and 7b resp. The breakdown voltage is 
defined here as the voltage at which a spark between the electrodes 
is observed if the voltage is slowly increased from zero to its final 
value. 

The spread in the observed breakdown voltages is influenced by 
the number of preceding breakdowns and by the time in which the 
voltage is adjusted from zero to the value at which breakdown 
occurs. 

Although the applied induction heating was sufficient to avoid 
discharges in gas liberated from the electrode surfaces, the main 
conditioning in order to obtain maximum breakdown voltages was 
done by previous discharges (see also § 4). These maximum values 
showed an average spread of about 8% for several measuring 
sequences. Even with well-conditioned electrodes, however, devi- 
ations to values as low as 40% of the maximum breakdown strength 
occasionally occurred. 

Electrodes not fully conditioned by previous discharges in general 
showed a lower breakdown voltage with faster increase of the 
applied voltage. Here also the pre-breakdown currents appear to 
_ have a conditioning influence. 


5.2. Pre-breakdown current. Contrary to the experiments 
described in § 4.2, the voltage in the a.c. experiments could be 
maintained during a longer time at an arbitrary value below the 
breakdown voltage. This facilitated the measurement of current 
and voltage, but on the other hand introduced phenomena which 
- influenced the pure vacuum breakdown, e.g., strong heating of the 
electrodes. 

Around the voltage where field emission becomes noticeable by 
the occurrence of peaks on the current records (fig. 10a), small red 
to white glowing spots appeared on the electrode surfaces. On 
raising the voltage, the number of these spots increased, whilst a 
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pale blue disk with a diameter of a few mm developed at the middle 
of each electrode surface. In the voltage interval just prior to 
breakdown these disks grew out to a transparent luminous channel 
connecting the electrodes. In this channel occasionally red glowing 
particles could be observed, crossing the gap. 

In most cases field emission from one electrode started at a lower 
voltage than from the other (fig. 10a), whilst initially the luminous 
phenomena also appeared on one electrode only. From the sign of 
the peaks on the current records it could be concluded that this 
occurred on the electrode emitting least. Apparently the glowing 
spots were caused by a bombardment of field-emission electrons. 
Some d.c. experiments confirmed this assumption. It should be 
noted that a strong X-ray radiation was measured during these 
experiments. 

At a much higher current (about 4 mA) some extreme cases have 
been observed, at which the whole of one of the electrodes was 
heated to a bright red glow without breakdown occurring. Exchange 
of electrodes proved this to be due to the “‘cold”’ electrode, which on 
the oscillograms also showed the largest field-emission currents. 


TPYTRY ee, 


a 


| ma [ray 
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Pig. 10. Pre-breakdown current and gap voltage for 50 Hz voltage. Copper 
electrodes, d = 0.1 mm. 

a) Peaks on the curve of current vs time caused by field emission (arrows). 

b) Current rising independently of the voltage. 


With a large series resistor (1.56 MQ) it was found that from 
certain values upwards the current could rise with decreasing 
voltages as shown in fig. 10b. It may be noticed that the current 
intensity through the gap measured in this case is larger than was 
found in the impulse voltage experiments just before breakdown 
(fig. 8). 


THE ELECTRICAL BREAKDOWN IN VACUUM 17 


This can be explained by assuming that some gas multiplication 
occurs in the vapour of an electrode hotspot, which, however, is not 
enough to lead to breakdown, firstly because the voltage is too low 
and secondly because there is no reason to suppose that in this case 
the current density will be high. This assumption is confirmed by 
the fact that the current may grow with decreasing voltage. With 
alternating current, hotspots on any of the electrodes can be formed 
during the voltage rise with less energy than with impulse voltage 
because the electrodes can be pre-heated by the preceding voltage 
cycles. 

This section so far has only dealt with experiments in the de- 
mountable vacuum chamber. In the all-glass system some com- 
parative experiments concerning pre-breakdown currents have 
been carried out. The same phenomena occurred here, except for 
some quantitative differences: 

a) At a given voltage the spread in the current values is larger 
in the all-glass system, 

b) At a given voltage the mean value of the current is larger in 
the all-glass system. 

It proved that the spread in the current-voltage curves is mainly 
due to taking data from experiments with different peak voltage 
values. A greater peak value giving a current-voltage curve that has 
shifted to lower values of the current. This is probably due to the 
more effective conditioning of the electrodes with high voltages and 
currents in the preceding voltage loops. Nevertheless, even if this 
effect is taken into account, the differences mentioned under a) 
and 8) still exist. The principal difference between the two set-ups 
might be the diffusion pump used. It is thought that oil-vapours 
resulting from the oil diffusion pump contaminate the electrode 
surfaces in sucha way as to lower the average effective work function. 

Finally the influence of the time interval after the induction 
heating of the electrodes has been investigated. Current data were 
compared from: 

a) cold electrodes, heated more than 24 hours ago. 
b) cold electrodes, heated less than half an hour ago. 
c) electrodes with a temperature of about 700°C. 

All data were taken with the same peak value of the alternating 
voltage. In the order of the different heat treatments mentioned 
above, the current intensity slightly diminished (at the same value 
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of the voltage). However, the differences obtained were smaller 
than the spread caused by taking values of the current at different 
peak values of the alternating voltage. 


§ 6. Discussion. Some conclusions deduced from the results dealt 
with in §§ 4 and 5 have been given in § 4.3 and § 5.2 respectively. 
Although our investigations cannot provide all information needed 
for the complete understanding of the mechanism of vacuum break- 
down in our experimental circumstances, it seems most probable 
that heating of the anode by field-emission electrons followed by 
evaporation of parts of the anode surface is the main factor initi- 
ating breakdown (§ 1, assumption @). 

It was observed that a luminous phenomenon and a high temper- 
ature at the anode is reached at a lower voltage than at the cathode. 
This is confirmed by observations with an ultra high-speed camera 
by Chiles 19), who also recorded the speed of the vapour stream 
leaving the anode (5 x 105 cm/s, impulse voltage experiments). 

It is less probable that heating of the cathode to a temperature 
where evaporation occurs should play a part in the initiation of 
breakdown (§ 1, assumption 6 and c). In that case one would 
expect that the projections which provide the field-emission 
electrons should always come into the same emission conditions, 
because these points are heated most. Consequently the different 
current-voltage relationships, given in fig. 7 should all coincide. 
However, it is repeatedly proved by our experiments that the 
effective work-function of the cathode tends to lower values for 
shorter rise times (shift of the curves in fig. 7 to the left). A 
possible explanation has been given in § 4.3 (removal of surface 
layers). 

From the observed luminous phenomena in a.c. experiments 
(§ 5.2) and from fig. 10d it appears that excitation and ionisation of 
vapour liberated at the anode already exists before the actual 
(sudden) breakdown takes place. The latter, then, might occur only 
when the vapour density passes a critical value. Fig. 10b might be 
explained by assuming that as the temperature of the electrodes is 
still increasing during the fall of the voltage, the vapour density 
also increases, thereby giving the current a chance to grow to a 
value limited by the series resistor. At higher voltages the vapour 
density might reach a value sufficiently high to start a multipli- 
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cation process resulting in a sudden increase of the current (break- 
down). 

The hypothesis that a critical energy flux to the anode is needed 
to initiate a discharge (§ 2.4) seems to be contradicted by the fact 
that, for electrodes not completely conditioned, breakdown occurs 
at low voltage and also at low pre-breakdown current (§ 4.2, table 
III). From a comparison of the current-voltage relationship in the 
a.c. experiments with the Fowler-Nordheim equation for the 
field-emission current density it appears, however, that clean 
cathode surfaces with a high effective work-function have a larger 
emitting surface than those with lower effective work-function. This 
means that on the less conditioned electrode surfaces there are only 
a few points that contribute to the emission of field electrons, thus 
giving a high current density, whilst the total current may be 
relatively low. 

As for Cranberg’s theory (§ 1, item d@), it seems unlikely that 
the observed irregular appearance of pieces of material alone is 
sufficient to cause breakdown at the well reproduceable voltage 
values found here. On the other hand it is possible that this theory 
holds for the casual occurrence of breakdown in large gaps, where 
at high voltages discharges can be initiated with irregular intervals 
of up to several hours. 
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A CLASS OF BOUNDARY VALUE PROBLEMS 
by W. E. WILLIAMS 


Department of Applied Mathematics, Liverpool University, England 


Summary 


A simple method is presented for the solution of the partial differential 
equation of diffusion type with constant values of the solution or its normal 
derivative prescribed on the surfaces of a wedge of arbitrary angle. It is 
shown that this solution may be transformed in such a manner that it yields 
the solution to a similar class of boundary value problems for the time- 
harmonic wave (i.e. Helmholtz’s) equation. Direct solutions are also 
obtained for the problem of the diffraction of acoustic or electromagnetic 
plane waves by a perfectly absorbing or a perfectly reflecting wedge. In one 
solution the difraction problem is solved by modifying the solution of a 
similar boundary value problem for the diffusion equation. In the second 
formulation the solution of the diffraction problem is obtained by expressing 
the total solution as a sum of a diffracted field and geometrical optics terms. 
The diffracted field is then obtained by imposing the conditions of continuity 
across the shadow lines of geometrical optics. The solution for the diffraction 
of an arbitrary plane pulse by a wedge is also obtained; the solution being 
valid even if the boundary conditions are of the impedance type. 


§ 1. Introduction. For over fifty years the solution of the dif- 
fusion and time harmonic wave (i.e. Helmholtz’s) equations with 
conditions specified on wedge surfaces has attracted the interest of 
a large number of investigators. Boundary value problems of this 
nature occur in a number of different physical situations, and many 
different methods of solution. have been considered. The object of 
the present work is to present a simple unified approach to the 
solution of a class of problems of this type. Before considering the 
techniques and problems examined in the present paper we shall 
review briefly the main methods used in the solution of boundary 
value problems on wedge and the physical situations which pro- 
duce such problems. 
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The problems for Helmholtz’s equation are mainly in acoustic 
and electromagnetic diffraction, the main interest having been 
concentrated on diffraction of plane waves by perfectly reflecting 
or perfectly absorbing wedges. The earliest solution is that of 
Sommerfeld 4) who employed the method of images on a Riemann 
sheet for a wedge of angle 2x. This approach was employed by 
Carslaw 2) who also obtained the Green’s functions for arbitrary 
wedges. An alternative approach to this type of diffraction problem 
is that of Macdonald 3) who obtained the solution as an infinite 
series of Bessel functions. Another different method is that of 
Kontorovich and Lebedev“) who employed integral trans- 
forms of Hankel and Bessel functions with respect to their order. 
This approach has also been utilised by Oberhettinger®) in 
deriving the Green’s functions for wedges. Diffraction problems 
for wedges with impedance conditions prescribed on the surfaces 
have been considered by the author §)’), the solution being effected 
by employing a suitable contour integral representation. This type 
of problem has also been solved by Senior 8) and by Maliuzhi- 
nets %) (the author has not seen this thesis). 

Another class of problems for Helmholtz’s equation which has 
attracted some interest is that when the solution or its derivative is 
prescribed on the wedge surfaces. The simplest case is when the 
prescribed values on the surface are constant; clearly one method 
of obtaining a solution is by employing the Green’s function. This 
latter function is rather complicated, and it seems desirable to use 
a more direct approach. Such an approach has been developed by 
Maliuzhinets 1%) who considers acoustic radiation from a wedge 
whose surfaces are oscillating with a constant velocity. His approach 
consists of using a contour integral representation and reducing the 
problem to solving an inhomogeneous difference equation. This 
latter equation is similar to one occurring in the problem of dif- 
fraction by a metallic wedge 6), and the solution is rather compli- 
cated. By employing certain properties of this solution Maliuzhi- 
nets shows that the solution for the radiation problem may be 
expressed as a real integral involving trigonometrical and hyper- 
bolic functions. The solution of the difference equation does not 
occur explicitly, and it thus seems that the introduction of this 
function is superfluous and that a more direct approach is feasible. 

The main boundary value problem for the diffusion equation 
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with constant values prescribed on wedge surfaces is the Rayleigh 

problem for a wedge in an incompressible fluid. Basically this 
problem is time dependent, but a Laplace transform reduces the 
basic equation to that of diffusion. A solution of this problem has 
been given by Hasimoto!!); his method involves a change of 
variable and a solution then obtained, by separation of variables, 
in terms of an infinite series of confluent hypergeometric functions. 
A suitable contour integral representation of these functions 
reduces the solution to a more suitable form. For wedge angles an 
integral sub-multiple of a a solution has been given by Sower- 
by 12), who employed the known Green’s function 2) for the prob- 
lem. This approach was subsequently generalized to a wedge of 
arbitrary angle by Sowerby and Cooke }8). The same type of 
boundary value problem occurs in the theory of heat conduction, 
and in this case a solution has been given by Jaeger 1}4) as an 
infinite series of Bessel functions. A similar boundary value problem 
occurs in the theory of neutron absorption, and a solution was 
obtained by Levine!) by means of the Kontorovich-Lebedev 
transform. A solution to this problem employing the Green’s 
function has been given by Lauwerier 16). 

In most of the above examples the final solution is relatively 
simple, though the intermediate calculations may be lengthy and 
involve complicated functions. It is thus clearly desirable to obtain 
a simple method of formulation which enables the above and more 
complicated problems to be solved in a simple and direct manner. 

In § 2 we consider the solution of the diffusion equation with 
arbitrary constant values specified on a wedge. The solution is 
effected by employing a real integral representation for the so- 
lution which was first used by Van Dantzig 1”). A simple appli- 
cation of the theory of distributions 18) or of generalised functions 19) 
reduces the boundary value problem to the solution of a difference 
equation. This latter equation is solved by means of Fourier 
transforms, the final solution being expressed in terms of trigono- 
metrical functions. The corresponding solution with constant values 
of the normal derivative prescribed on the wedge surface is also 
given. The approach may be applied to the case when a linear 
combination of the solution and its normal derivative has pre- 
scribed constant values on the two wedge surfaces. It may also be 
used to solve problems where such a linear combination vanishes 
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on one wedge surface and the function is constant on the other 
surface. Problems of this type occur in oceanography 2°). It is also 
shown that the solution obtained may, by slight manipulation, be 
written in such a manner that it may also be applied to obtain 
solutions of Helmholtz’s equation. Thus solutions are immedi- 
ately obtained to the type of problem considered by Maliuzhinets, 
and the solution obtained by the present approach agrees with that 
of Maliuzhinets. 

In §3 the particular type of boundary value problem occurring 
in acoustic and electromagnetic diffraction theory is considered. 
Two different methods of solution are examined. The first method 
is an indirect one; the approach of § 2 is used to obtain the solution 
of a particular type of boundary value problem for the diffusion 
equation. It is then shown that, by suitable re-arrangement and 
slight modification, the solution for the diffraction problem may be 
obtained. In the second approach the solution to the diffraction 
problem is formulated as a sum of the plane wave terms of geo- 
metrical optics and a diffracted term. A real integral representation 
of this term is employed, and the conditions of continuity across 
the shadow lines of geometrical optics yield a set of simultaneous 
functional equations. The solution of these is obtained as in § 2. 
The approach may also be used to solve problems where impedance 
conditions are prescribed on the wedge surfaces. The only difference 
in these problems is that the functional equation to be solved is 
more complicated. 

The particular form of the solution obtained by this second 
approach is particularly suitable for the solution of diffraction by 
plane pulses of arbitrary shape. The solution for the time-dependent 
problem may be written down immediately from the time-harmonic 
solution if the boundary conditions on the wedge surface are of the 
form (time derivative on surface/normal derivative on surface) = 
constant. This type of problem has been considered by Papa- 
dopoulos #1), for the particular case of a half-plane, by assuming 
dynamic similarity. It is also of interest to note that the diffraction 
problem for the time-dependent wave equation may be formulated 
directly without recourse to the time-harmonic solution. The so- 
lution is once more expressed as the sum of geometrical optics 
terms and a diffracted term. The only difference lies in the integral 
representation of the diffracted term; the formal steps in the 
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_ analysis and the difference equations are exactly as for the time 
harmonic case. This approach seems to present a more direct and 
simpler approach to these problems than that of Keller 22) or 
Papadopoulos 23), It also has the advantage that it is not re- 
stricted to unit step function pulses. 


§ 2. Simple solutions of the diffusion equation. Throughout the 
present work we shall be concerned with the region of physical 
space outside a wedge of angle «. If (7, #) denote cylindrical polar 
coordinates with axis parallel to the edge of the wedge, the region 
of space of interest is defined by. O= 7. = oo, 0 = 0 =< ce. In the 
present section we shall obtain the solution of two simple boundary 
value problems for the equation 


(V2 — k®)p = 0, (1) 
where 

e oped Rs Bees 
y OF 25 v2 of? 


and k is real. The first problem is the determination of a solution 
mi. wey Or (1) such that ¢) = 1 one = 0, di = 0 on. 0 = a and 
$1 >0 as y +>oo in 0 <#<a. The second problem is that of 
obtaining a solution ¢o(7, 7) with 7v—lé¢2/ee = 1 on 8 =O and 
71é¢2/00 = 0 on & = «; ¢d2 satisfies the same conditions as ¢j at 
infinity. Clearly from ¢1, 2 we can construct a solution of (1) with 
values cj, co on # =O and « respectively, such a solution is 
cipilr, B®) + cobi(7, « — #). From these functions ¢1, ¢2 solutions 
can be obtained for the Rayleigh problem for an incompressible 
fluid 11) and for problems in heat conduction treated by Jaeger 14). 
In both these problems # will involve a Laplace transform para- 
meter. The solution to the diffusion problem treated by Levine }) 
may also be obtained. We shall not consider the application of our 
results to these particular problems but merely obtain the solutions. 
These solutions will be obtained by employing an integral 
| representation for solutions of (1). The most convenient for the 
present problem is that of Van Dantzig 1’). We thus have 


$1 = fete {Fily + i(8 — «)] — Fily — 18 — a}]}dy. (2) 


F(n) is an arbitrary function, analytic in the strip —« < Imy <~«, 
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and the form of ¢1 is such that the conditions on 0 = & are auto- 
matically satisfied. The condition on # = 0 gives 


fe at (1p tee) fae oN fae (3) 
The part of the integrand in square brackets in (3) clearly does 
not represent a function in the normal sense, but nevertheless it 
may be interpreted as a generalised function 1%) or as defining a 
distribution 18). In either case its Fourier transform may be defined. 
For convenience the symbolic notation of the delta function will be 
employed and thus 


Py(y — to) — Fi(y + ta) = d(q). (4) 


Equation (4) may be solved by conventional Fourier transform 
techniques. If F'1(¢) is defined by 


then (4) becomes 


2F,(C) sinh af = 1. 


Thus 
1 by } - sin 7 
IB aG =5— [or \Ne= | 1 ha el 
1() Do le 1(¢)dé che d¢ = b+ tanh der, 
where 7 == /¢. Hence 
t kr sinh », . 
pi = me e~*rsinha Toth 47(y + 78) — coth 47(y — i8)]dy. (5) 

Similarly 


a ikr sinh » } 
bo =z e7mran®" [coth $r(n + 79) + coth 47(y — i8)]dy. (6) 


—oo 


We now apply these results to one of the boundary value problems 
mentioned above, namely that treated by Levine. In this case the 
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function ¢ is required to be unity on both surfaces and hence 


C= an e *rsinhy (cothtr(n + 18) — coth $7(y — 1) + 
+ coth $7(y + t« — 7) — coth $r(n — ia + i9)]dy = 


es e7*rsinb"’cosech r(n + 10) — cosech r(y — i8)]dy. (7) 
lod 


Eq. (7) is identical with the result obtained by Levine; also with 
a suitable interpretation of the parameter k it represents the solution 
for the Rayleigh problem for an incompressible fluid and a wedge of 
exterior angle «. It should be noted that (7) could have been 
derived more directly by employing the properties of symmetry 
which require that 64/20 = 0 on # = 3a. 

It is clearly of interest to consider whether the above results may 
be employed to obtain solutions of 

Ve) oO (8) 

with either ¢ or its derivative having prescribed constant values on 
the wedge surface. Boundary value problems of this nature occur 
in acoustics and electromagnetic theory in determining solutions of 
the time-harmonic wave equation with prescribed conditions on a 
boundary. Formally at least one would expect to obtain a solution 
of (8) from that of (1) by replacing k by tk. Clearly equations (5) 
and (6) are such that the integrals diverge for k purely imaginary, 
and thus it is necessary to re-write these in a form which is valid for 
such values of k. 

We have that 


cot bani 
4 


br =~ | e~P-*[coth dr (—i4-+ 18) —coth}r(y— dia —i8)]dy— 
(04 
“ —oo+ tani 
cot knit 
=) | eem(coth dr(y—ah+ 10) coth 4r(n—din—18)]dn — 
lo 4 
—oco+ $i 
co— tri 


Pe | e~*reoshn(eothd7(n--ida—id) —coth$r(y+ 310+ 19) | an} . (9) 


—oo— $n 
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Both contours of integration in (9) may be deformed to the real 
axis; some poles of the integrand may be captured in this deforma- 
tion. Thus (9) becomes 


[5(—9)]7 [Laat o]- 
p1 a S exp[— kr sin(# + 2nx)] — Y exp[— er sin(@ — 2na)|-+- 
n=0 peas 
=e oe jPeres sinh Ay - vd) dy. (10) 
2a cosh 7(7 + 73) — cos 4ar 


0 
The notation [x]+ refers to the nearest integer less than, or greater 
than, x respectively. 

Similarly it can be shown that 


i [5-G—-0)) 


o2 = Al > exp[— kr sin(? + 2na)] + 


=0 


[sGa+a)]- 
+ ¥ exp[—Arsin(d — os = 
n=[=]* 


e —kreoshn 


— sin d : 
Se | cosh r(7 + 70) — cos $ur 
0 


dy. (11) 


k may now be replaced by 7k in (10) and (11), and the functions so 
obtained wiil represent solutions of (8) satisfying the same boundary 
conditions as ¢; and ¢z on the wedge surface. The resulting func- 
tions will, apart from the plane wave terms, satisfy the radiation 
condition 
lim 7?(@6/ér + ik) 

and thus, in terms of solutions of the time-harmonic wave equation, 
will represent outgoing waves with time dependence exp dkct 
where c is the wave velocity. The solution appropriate to the time 
dependence exp(— tkct) is obtained by replacing k in (10) and (11) 
by — 7k. 

Maliuzhinets 1) has considered the solution ¢ of (8) with 
v-10¢/69 = 1 on &=0 and 06/08 =0 on & =a and, when no- 
tational changes are taken into account, the function obtained by 
replacing k by tk in (11) agrees with Maliuzhinets’ solution. 
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The method used is directly applicable to problems where im- 
pedance conditions are prescribed on the wedge surface, the only 
difference being in the form of the difference equation (4). For 


example if we have ¢; = 1 on 8 = O and an impedance condition 
prescribed on # = «, (4) becomes 


Fi (y — ta) — G(n)Fi(n + tx) = d(n). (12) 


G(y) is a known function of 7 and the impedance. The solution of 
(12) is 


tH (x 
Fy(n) = Wis tanh 477, 

where H(7) is an appropriate solution of (12) with the righthand 
side set equal to zero. Solutions of this equation appropriate to 
impedance conditions occurring in acoustics and electromagnetic 
theory have been considered by the author §)’) ; a similar equation 
has been encountered by Peters 24) and Van Dantzig 1%). 

In the following section we shall consider whether the techniques 
of the present section are applicable to diffraction problems. 


§ 3. Diffraction by a wedge. The basic problem in diffraction 
theory is the determination of a solution of (8) satisfying certain 
conditions on the wedge surfaces and at infinity. The condition at 
infinity depends on the incident field; in the present case we shall 
consider the incident field to be that of a plane wave. The condition 
at infinity is that the total solution, apart from the terms predicted 
from geometrical optics, satisfies a radiation condition. We shall 
employ the radiation condition appropriate to a harmonic time 
dependence exp zkct. The conditions at the wedge depend on the 
physical problem considered. For example if ¢ represents the 
acoustic velocity potential, then 44/20 = O represents the condition 
on a perfectly hard surface and ¢ = 0 the condition for a perfectly 
soft surface. If an impedance condition is to be prescribed, then 
rp/(é¢/@9) is constant on the wedge surface. Since the object of the 
present work is to indicate a simpler and more direct solution of 
these problems, it is sufficient to consider only ¢ = 0 on the wedge 
surface. The solution for electromagnetic diffraction with im- 
pedance boundary conditions on the wedge surface has been ob- 
tained by a slightly different approach by the author 6\7\,. 
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The incident wave is taken to be exp kr cos(# — 8) and the 
total solution 4 written as exp ikr cos(® — Bo) + 7. One of the 
difficulties of diffraction theory is that it is not possible to obtain a 
real integral representation for z. It is, however, possible to formu- 
late the diffraction problem rigorously by writing ¢ as the sum of 
geometrical optics terms and a diffracted term which may be 
represented by a real integral and which is discontinuous across the 
shadow lines of geometrical optics. It has in fact been pointed out 
by Peters and Stoker 25) that diffraction problems should be 
formulated in this way in order to obtain unique solutions. We shall 
eventually consider the solution of the diffraction problem from 
this rigorous viewpoint. 

An alternative approach which we shall now consider is de- 
termining the solution of the diffraction problem from an appro- 
priately formulated boundary value problem for (1). The advantage 
of employing (1) is that it is often possible to obtain a real integral 
representation of the form of (2). We therefore consider a solution 
yy of (1) with yy = — éxp(—2e7 sinh 44) on =O oe 
= — exp(— tkrsinh dg) on & = «, Ay, Ag real. The approach of 
§ 2 gives 

a I —ikrsinh 1 c 4 
A e [coth 47(n + 79 — a1) — tanh 47(y + 76 — Ag) — 


—co 


— coth $7(7 — 7 — dy) + tanh $7(y — 19 — Ag)]dy. (13) 


Clearly from a formal viewpoint one would expect to determine the 
appropriate solution x of (8) from (13) by replacing k by zk and 
Ay, Az by (3 — Yo), 1($ + Bo — «) respectively. We consider first 
the form of (13) for 41, Az having these values; for simplicity the 
following discussion will be restricted to the case of « of z+ Bo. 
Any other value of « may be similarly treated. 

For the above imaginary values of 4; and 4» it is clear that val 
vanishes on # = 0 and « and also the integrand is singular when 
& = 4a + 0, and we therefore consider the behaviour of yintieat, 
these critical values of #. Clearly this reduces to considering 


1(6) = ferirrinby coth 37(n — 10)dy 


—oo 


for small 6. If the infinite series for coth z is substituted into the 
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expression for /(6), then we obtain 


Ge eT 20) eet Art 


Thus, as # increases from 0 to «, y; decreases by | across the line 
3 = 42 — 9 and increases by 1 across # = ix + Hp. Thus a 
continuous solution yz may be obtained from y1 by defining x2 by 
41 — exp kr cos(# + Vo) 0<%< tna — %, 
2 = 4971 an —to << An + Ho, (14) 
| x1 — exp kr cos(8 — Jo) = 4x +H < 2B. 
Furthermore 


72 = — exp(Arcos Jo) on d= 0 and v2 =— exp[krcos(a—)] on P=. 


It thus appears that a solution to the diffraction problem could be 
obtained if v7; were re-written in a form valid for purely imaginary . 
From (13) and (14) by a method similar to that used in deriving 
(10) and (11) we obtain 


Fs =o = 00) SF (+ Ho — 8) + 
rt OD too) — FY, no — 0) + G(r, 8), (15) 


where 


em kr cosh 


dy (16) 


] 
Fir, 6) = —sin 76 
2a 


0 


cosh 77 — cos 70 


and 
( — exp kr cos(# + Io) O<9<a— Vo, 
Gy, 2) =} 0) z— 99 <P <a+%, (17) 
| — exp &r cos(8 — Ho) F>x+ Do. 


The solution to the diffraction problem may now be obtained from 
(15)—(17) by replacing & by 7k, and the solution so obtained agrees 
, with that of Oberhettinger ”§). 

A rigorous formulation of the diffraction problem will now be 
examined. The solution ¢ may be written 


exp ikr cos(9—%o) —exp tkr cos(9+-9o) +f1 0<b<a—Vo, 
exp tkr cos(# — Jo) + fe ar—09<9<2+, (18) 
fs O>24+-Vo0. = 
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The functions /1, fz, /3 are required to satisfy the radiation condition. 
Clearly appropiate integral representations for f1, fa, fs which 
satisfy the boundary conditions on the wedge surfaces are 


f= fexp(— thr cosh n)LFiln + #) — Fin — Olde 


fem fexpl=dhF cosh LEG eN Geece 
fg = fexp — thr cosh yi Fa(y 4 40 — ta). Page tor tee 


The conventional methods of asymptotic expansions 2”) show that 
the above expressions will satisfy the radiation condition if the terms 
in square brackets are absolutely integrable. The function ¢@ is 
required to be continuous everywhere and hence 


Pily -- 10 = 9) = Fil — 4 0) — 8) 
= Fly + 1(% — 9o)] — Galn — 1(% — Fo)], 
F ily + .1(% — Oo)] Pity — te — Yo) | = 
= Fe[n + ua — Bo)] + Galn — 1% — Do], 
Fely + t(@ + 8o)] — Galy — + + Bo)] + 4(m) = 
= F3[n + t(n + Bo — a)] — Fain — 1(a + 80 — @)], 
Faln + t(% + 8o)| + Galn — 1 + Bo)] = 
= Faly + t(a + Bo — «)] + Fly — t( + Bo — a)]. 


The above simultaneous difference equations are soluble by trans- 
form methods as in § 2, and it is easily found that 


— 4tFi(n) = coth $7[y + ia — 189] — coth 47[m + i(a + B)] + 
+ coth $7[y + 1(8 — x)} — coth 47[7 — i(a + Bo)]. 


The solution obtained agrees with that derived from (15)—(17). 
The above method of formulation may be applied to the boundary 
condition 0¢/09 = 0 and also to problems with impedance con- 
ditions prescribed on the wedge surfaces. The solution in all these 
cases splits up into the terms of geometrical optics and a diffracted 
term which has an integral representation similar to that of hivtawiss 
In particular if the impedance on the wedge boundaries is such that 
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on these boundaries @/00/¢ = ikZr where Z is independent of k, 
then the integrand for the diffracted field only contains k as 
exp — tkr cosh 7. It will now be shown that this fact together with 
the form of the solution obtained by the rigorous formulation 
enables the solutions for diffraction by an arbitrary time-dependent 
plane wave to be written down immediately. For impedances of 
the above form we have that, for an incident plane wave, 
exp 7k{ct + r cos(# — %)] the solution ¢ is of the form 


Plane wave terms im 
e'*et 1 of geometrical + (Sa F(n, 0)dy $, 
optics 0 


where ¢F/@k = 0. The corresponding solution for a plane pulse 
Gict + r cos( — Ho)] with 4/00 = (Zr/c)(éd/et) on the wedge 
boundaries is thus given by 


co 


Sum of plane 
pulse + | G(ct — rcoshy)F(n, &)dn. {19) 
terms 0 
Boundary value problems of this type have been investigated by 
Papadopoulos 2!) using arguments of dynamic similarity. The 
present approach seems to be far more direct and also gives so- 
lutions for arbitrary incident pulses and not for those of step 
function form only. The quantity evaluated dy Papadopoulos 
is é¢/ét, and for a unit step function the contribution to this from 
the integral of (19) is 
F(cosh-1 ct/r, 0)H (ct — 7) 
where H is Heaviside’s unit function. The contribution is thus 
limited to the region ct > 7 which corresponds to the region in 
which the wave equation, re-written in terms of the variables 
(r/t) and 9@, is elliptic. For ¢ or @$/@3 zero on the wedges the so- 
lution (19) has been given by Friedlander *8). 

It is of interest to note that these solutions of the wave equation 
could have been obtained directly without recourse to Fourier 
superposition. The time-dependent diffraction problem could have 
been formulated as in (19) with the functions /1, fe, fg having re- 
presentations of the form 


fGlet —rcosh y)[Fi(y + 18) + Fa(n — vd)dn. 


—co 
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The whole analysis would then proceed as before. This approach 
seems to yield a more general method to the solution of time- 
dependent diffraction problems than arguments employing dynamic 
similarity or the conical flows method of Keller and Blank 22). 
The method seems much simpler in that a direct solution of the 
difference equations is obtained, and it is unnecessary to ‘spot’ 
analytic functions possessing appropriate singularities. 


Received 14th July, 1960. 
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ESTIMATING THE ELECTRIC FIELD INSIDE A 
RECTANGULAR TANK WITH BOUNDARIES 


AT ZERO POTENTIAL 
by S. J. VELLENGA 


Bataafse Internationale Petroleum Maatschappij, the Hague, Netherlands 


Summary 


The maximum field strength occurring in a tank of rectangular shape 
whose boundaries are at zero potential and which is partially filled with a 
liquid bearing a constant charge density can be quickly calculated by means 
of a strongly convergent series derived in the article. Restrictive assumptions 
are: 1. The charge is homogeneously distributed, 2. The difference in the 
values of the dielectric constant of the vapour and the liquid is ignored, 
3. The surface of the liquid is parallel to the bottom of the tank. 


List of symbols 


ab. 
t 


= 
S 
3S 


3 


MRaAMHAE TSS 


length, width and height of the tank, respectively (m) 
electric current entering the tank (4) 

summation indices 

B2m? + n2 

sV = total charge in the liquid (C) 

charge density (C/m?) 

time 

fraction of the tank’s height filled 

coefficient in the series for E 

dielectric displacement (C/m?) 

field strength (V/m) 

potential (V) 

c/a 
c/b 
fraction of the tank’s height at which the field is measured 
relative dielectric constant 

absolute dielectric constant 8.854 x 10-12 As/Vm 

a/b = b/c; shape parameter 

relaxation time (s) 


shape parameters 
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§ 1. Introduction. The problem discussed in this article was 
prompted by the author’s interest in the phenomenon of static 
electricity as it occurs in the oil industry. The underlying equation 
(Poisson’s equation) occurs in other fields of science and therefore 
the derivation given below may be of interest to others as well. 

Pumping a hydrocarbon liquid through a pipe or a filter often 
causes the liquid to become electrically charged. If such a charged 
liquid passes into a tank, the electric field associated with the 
charge will become noticeable and may increase sufficiently in 
strength to cause an electric spark discharge. As this presents a fire 
risk, it is desirable to be able to predict the magnitude and distri- 
bution of such an electric field. 

For a tank of given shape the electric field is determined by the 
amount and distribution of charge in the liquid and by the size of 
the tank, but the shape of the tank also greatly influences the field 
strength. The present article discusses the effect of the dimensions 
of a partially filled rectangular tank on the electric field inside such 
a tank. The simplifying assumptions are: 

a. the boundaries of the tank are at a constant potential, 

b. the charge is homogeneously distributed, 

c. the difference in dielectric constant between vapour and liquid 
can be ignored, 

d. the surface of the liquid is parallel to the tank’s bottom. 

In practice nearly all tanks are made of metal, and condition a is 
therefore usually fulfilled. A homogeneous distribution of charge, 
while less common, may exist in tanks which are being filled at 
high speed, and fortunately this is a case of great practical interest 
from the point of view of safety. Condition c is of course never 
realized in practice, but since the dielectric constants of vapour and 
a hydrocarbon liquid differ only by a factor of about 2, the approxi- 
mation implied by condition ¢ is not very severe. 

Of course, if the tank is nearly full, the field strength in the liquid 
will no longer be affected very much by the dielectric constant of 
the vapour. Since the field strength remains finite, the surface of 
the liquid when about to touch the roof will be equipotential with 
roof and walls. Consequently the field in the liquid is not influenced 
by the presence of the thin air gap at the top. In the absence of 
surface charges the dielectric displacement is continuous across the 
liquid surface so that the field strength in the air gap is proportional 
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to the field strength in the liquid at the surface, in inverse ratio to 
the dielectric constants. 

The statement of the problem in the above form was first given 
by W. E. Matheson and co-workers in an internal report of the 
Douglas Aircraft Company Inc. (see 1)). On the basis of their work 
the present article gives a rapid method for the calculation of the 
maximum value of the field strength in a tank as described above. 


§ 2. Tanks of large lateral dimensions. The electrical situation 
in a tank whose height is small in comparison with its lateral 
dimensions can be compared to that occurring between two in- 
finitely extended flat plates, separated by a layer of charged liquid 
and a layer of vapour. 

This case has been fully treated in reference 2), p. 104, where it is 
shown that allowance can be made for the difference in dielectric 
constant between liquid and vapour by means of a correction 
factor. As a first approximation this correction factor can of course 
be applied to the equation for the field strength derived below 
(equation (2)), particularly if the level of the liquid is very close to 
the top of the tank. 

In that case the situation in the tank far from the walls approaches 
that between two infinitely extended flat plates. Finally, it is also 
shown in 2) that there may be an additional charge on the surface 
of the liquid. However, as this surface charge is not so likely to 
develop in the well-stirred tanks considered here, it will be ignored. 

The relationship between the field strength E and the height c of 
the tank is the following: 


pes (1 


s = space charge density, 
w — fraction of the height to which the tank is filled, 
e = relative dielectric constant of the liquid and the vapour, 
€9 = absolute dielectric constant of vacuum. 
The notation differs slightly from that of ®). 


§ 3. General treatment. If the lateral dimensions of the tank are 
not very large with respect to its height, allowance must be made for 


38 S. J. VELLENGA 


them. This usually involves considerable computation, but for the 
special case considered here the arithmetic can be reduced ap- 
preciably. 

For a completely full rectangular tank, with boundaries at zero 
potential and a homogeneous charge distribution, the potential 
field found by the method of Green’s functions is the same as the 
temperature distribution in a rectangular solid with a uniform 
distribution of heat sources and the surface kept at zero temper- 
ature. This field is given in 3), page 1259. W. E. Matheson has 
extended the result to the case of partially filled tanks (see !)). 

Both 1) and 3) give the potential distribution in the form of a 
triple summation, each summation extending to infinity. This form 
is of course not very convenient for assessing a given situation. In 1) 
the difficulty is met by considering only the vertical field strength 
in the centre line of the tank, where it is larger than anywhere else. 
It is further assumed that the width of the tank is the geometric 
mean between its length and its height. By this means one shape 
parameter is removed from the complete equation. For a number of 
values of the remaining parameter the triple summation was 
performed over 729 terms, using an IBM 701 computer. It thus 
became possible to plot the field strength against a dimensionless 
number characterizing the liquid. 

The field strength in the vertical centre line of a tank can, how- 
ever, also be represented as a single series which converges rapidly, 
provided the height of the tank is not too much larger than its 
lateral dimensions. In this equation both shape parameters appear, 
the assumption regarding the relation between width, length and 
height of the tank being unnecessary. This equation reads as 
follows: 


SC me 
Fees} or w2 — ~ An(cos anw — 1)cos any |, (2) 
where 
Ay =a numerical factor depending on the shape parameters. 
Values for these factors are give in the appended table and 
graph. . 
y = the fraction of the tank’s height at which the field is 
measured (y > w). 
The derivation is given in Appendix I. 
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Fig. 1. Correction terms A, in the equation for the field strength in 
rectangular tanks. 


The main difficulty in applying (1) and (2) lies in finding a 
representative value for the charge density s. Ref. 2), p. 97, gives for 
a full tank 

ats =r, (3) 
where 

a, b, c = length, width and height of the tank, respectively, 

i = electric current entering the tank, 
+ = relaxation time of the liquid. 


Strictly speaking, this equation refers to an electrically stationary 
condition, but in practice this is approximated fairly well, because 
the time of filling is usually much larger than the relaxation time. 

Equation (3) can also be applied if the tank is not wholly filled. 
If necessary, a correction might be made according to the line of 
thought given in 1), p. 99. 

Before (2) became available, the opportunity presented itself of 
testing the general approach outlined in the above by means of 
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some experimental data. As it was not considered to be worth while 
setting up a computer programme for this single case, the necessary 
summation of the triple series was performed with the aid of a desk 
calculator, and took about three days. The difference between the 
measured and the calculated value was only 4%. In this particular 
case the tank was nearly full and a correction was made for the 
difference in dielectric constant. 

Subsequent recalculation of the field strength by means of (2) 
took only about half an hour; with the aid of the appended graph 
the calculation is only a matter of minutes. 


APPENDIX I 


Derivation of the field in a rectangular, earthed tank, containing a 
uniformly distributed space charge 


Green’s function for this case is given in 8), p. 1258. It should be 
noted that the numerical factor given in this reference is low by a 
factor of 4. Moreover, Green’s function as given should be divided 
by a factor of 4% ee9 if a rationalized system of electrical units is 
used. The coordinate system used has its origin at a corner of the 
box. The box is thus bounded by the x, y and z-planes and the 
corner opposite the origin has the coordinates (a, 6, c). Green’s 
function is 


8 Ss 
G(r/ro) = Lia)? = =) 
(0) = —— Ella) + (mb)? + (nfo)2} 
_ wy . aMYo . Neo . wx | amy . mMNz 
. sin Sie ees sin sin. —— sn ° 
Cc a b Cc 


where r is the vector from the origin to the point (x, y, z) and the 
subscript 0 refers to the location of the point source. 

In order to obtain the total potential distribution in the box the 
integration must be carried out over all charged points in the 
liquid volume. If the liquid level is at cw, and since we have as- 
sumed the charge density s to be constant throughout the liquid, 
only the following integration need be performed: 


ab p 


a . MX. EM 
sin sin -— 5 — sin ——— > dxodyodzo. 


00 0 
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Hence the potential distribution is given by 


= (32 a ( l i 4 ( ik 3 ( n We (1 — cos anw) 
go ““es5 1 a b NG lmn 
wx. . my | mMNzZ 


Sti sin sin 
a b é 


, 


where the summation indices / and m only take odd values, but ” 
takes all integer values. 


Vertical freld strength 
The vertical field strength Ez is found by differentation 
1D SSS ou /az. 

The maximum value of E£;, seen as a function of x and y, occurs 
at x = 4a, y = 4b, because there 0E,/0x = O and 0E;/dy = 0. This 
maximum value of £; is still a function of z, and in principle 
@E,/éz =O should be formed in determining the location of the 
absolute maximum of E;. However, it is evident physically that the 
maximum field strength £, in the vapour space will occur at the 
liquid’s surface. Since the latter value is the most interesting and 
can be readily computed from equations as given here, it is not 
necessary to write out 6Ez/éz = 0. The maximum value of Ez, asa 
function of z, is now found to be 


Se BO Oe, L\2 m \2 n \2]-1(1 — cos anw) 
; See Ss a aa eS aioe 
m4 eepC 1 a b c lm 
& ls et ANZ 
. sin —— sin —— cos ; 
2 2 c 


where it is no longer necessary to specify odd values for the sum- 
mation indices / and m because the even terms become zero in any 
case. This triple series can be summed with respect to each of its 
summation indices in turn because it is absolutely convergent; 
see 4) and 5). 

The summation over / can be written as 

sin _sinlg 
Ua? +p?) 
where «= c/a; p2 = Bm? + 77; B= c/b, whilst m will in due 
course assume the value $x. Now 


oe 


i 
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ie, 1 5 cos lp ty 
dp) a 1 P+ pat 
1 [ aa cosh(x — 9)p/« a ] 
Te E sinh wp/« 2p2 I 


For the derivation of this summation see 6). Integrating this 


result, 
eee! ee ( sinh(a — ¢)p/« ) aa 
~ a L 262 sinh xp/c 2p2 04 


where the constant of integration is found from the condition 
= 0 forig = 0. hog eZ 


I= rip E es | 


Substituting this result and writing y = 2/c, 


ee 2 We ee I E 1 ] 
w+ esq 1 4m(B2m2 + n?) cosh(pzt/2a) 
.(1 — cos mnw)sin am cos any. 
The double summation can be separated into two series, one of 
which is 
1A 


4m/(B2m2 + n?) 


x 
1 


(1 — cos anw)sin 42 cos any. 


In this series the summation with respect to m can be carried out in 
the same way as has been done for /: 


I sin mp 702 E l ] 
4 1 m(p2m + 2) 16n?2 cosh(nz/2B) 1 
Substituting this result, 


32 a 2 
eae SC S [ It (: of l ) i 
m4 ee 1 L 16n2 cosh (nz/26) 
z sin ms/2 
1 4mp? cosh pa/2« 


The last factors can be written as a sum of cosines: 


] (1 — cos anw)cos any. 


cos any(1 — cos anw) = 


= — [} cos an(y + w) + 4cos an(y — w) — cos any |. 
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TAB EBT 


Coefficients A » 


4 1 16} ~e2 sin (mz/2) 
A n= - + > = 
a nN 3G LT ; Daw 
n2 cosh | — — Mp cosh | ———— 
B 2 metre 
B a} 0.00 0.25 0.50 G75 1.00 1.25 1.50 fore) 
0.25 0.002 0.003 
0.50 0.035 0.035 0.059 | ne il 


0.75 0.098 0.098 0.111 0.146 
1.00 0.162 0.162 0.168 0.188 0.217 
1.25 0.213 0.213 0.216 0.228 0.246 0.266 
1.50 0.253 0.253 0.254 0.261 0.272 0.286 0.300 0.405 


BN a} 0.50 0.75 1.00 1.25 1,505) 1.60. 9] 


0.75 0.000 0.000 
1.00 0.001 0.001 0.002 
1.25 0.002 0.002 0.003 0.004 
1.50 0.004 0.004 0.004 0.005 0.007 0.045 


1.00 
1.00 0.000 
1.25 0.000 0.001 n=4 
1.50 0.001 | 0.001 0.001 0.002 


1.25 


1.25 
1.25 | 0.000 
1.50 | 0,000 | 0.000 
Now 
>> Staaed = }(02 — 210 + 3m?) (see e.g. ref. *), p. 3). 
i OM 
Therefore 
> uae cos zn(y + w) + cosan(y — w) —2cosany _ Lae ie 


1. 16 2n2 See A 
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Substituting this result, 


E Ss ie: 4 s (cos wnw — 1)cos amy _ 16 s Ss 
i n2 cosh nz/2B TE al 


(cos mnw — 1) cos any sin ses 
mp? cosh pa/2a 


The series in this expression usually converge extremely rapidly, 
so that in most cases a few terms suffice to give a reasonable 
approximation. Slight rearranging and numerical evaluation of 
certain constants lead to (2) in the main part of this report. The 
factors A» can be calculated from 


4 1 16 sin mz/2 
x2 n2 cosh nz/2p x3 ~ mp? cosh pa/2a 


A == 


This series converges very rapidly indeed. For example, if « = 0.724 
and 6 = 1.484, the first term is 0.2510 and the second term is 
already extremely small: 0.00661. By double linear interpolation in 
table I a value of 0.258 is found. 


Received 23rd July, 1960. 
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Summary 


The two-dimensional stagnation flow of an electrically conducting, 
incompressible and viscous fluid against a plane wall is investigated for the 
case when the induced field is negligible compared to the imposed transverse 
magnetic field. It is found that the component of the velocity parallel to the 
plate as well as the drag coefficient decrease with the increase in the magnetic 
field. Furthermore, it is observed that the velocity component parallel to the 
plate is essentially constant except in a layer of constant thickness, a result 
which is true in the non-magnetic case also. 


§ 1. Introduction. The problem of two-dimensional stagnation 
flow of an incompressible viscous fluid against a plane wall was 
completely solved by Blasius +) and Hiemenz 2). Blasius reduced 
the partial differential equations for the flow to ordinary differential 
equations with requisite boundary conditions, and these were 
integrated numerically by Hiemenz. The solution of Hiemenz 
was, however, improved by Howarth 8). This exact solution has 
close relation to solutions of the boundary-layer type for flows past 
an arbitrarily shaped body and substantiates the assumptions of 
the boundary layer theory in the most important region of the field 
of flow. The present note is an extension of this problem to the case 
of an electrically conducting fluid in the presence of a transverse 
magnetic field. This solution, however, has an important bearing 
on missile re-entry problems. When a missile re-enters the earth’s 
atmosphere, a sufficient amount of heat is generated due to friction 
of the air molecules. This viscous heating may sometimes be so 


re Se Te 
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large as to ionise the air near the forward stagnation point. As 
ionised air is electrically conducting, a magnetic field may be 
applied to induce forces in the air which in turn will be retarded, the 
rate of heat transfer at the wall being reduced as a consequence. 


§ 2. Basic equations and their solutions. The stagnation point in 
two-dimensional laminar flow against an infinite plate is taken as 
the origin, y is the distance from the plate and x the distance along 
the plate in the plane of the flow. A constant magnetic field Hp is 
applied along the y-direction. 

The equation of continuity is 


Leona lr toca (1) 
Ox oy 
The equations of momentum are 
o ou Yaa 02u o2u ob? 
flea Sea p > \= u, (2) 
Ox oy p OX Ox? oy? p 
0 dv Teg 02y 02u 
ee +v Ss P ae ( ) , (3) 
Ox oy p oy 0x2 oy? 


where o is the electrical conductivity and B = weHo, ue being the 
magnetic permeability. In writing these equations, the secondary 
effects of magnetic induction are ignored 4). We have also assumed 
that the electric field is zero; this seems to be a justifiable assump- 
tion since no external electric field is applied and the effect of 
polarisation of the ionized field may be expected to be small if two- 
dimensional conditions occur in the ionized layer 5). In view of 
these assumptions, the Maxwell equations become superfluous. 

The velocity distribution in the potential flow in the neighbour- 
hood of the stagnation point at « = y = 0 is given by 


Uco = BX; Vo = — By. (4) 
Using (4), the pressure gradient @f/@x is now obtained from (2) 
with y = 0 as 


1 @ 
— =m, (5) 
where 
B2 
p? = pe + 8. (6) 
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It follows from (5) and (6) that the pressure distribution is affected 
by the magnetic field. 
If the stream function yp is used such that 
ow oy 
4=—-, v=——, (7) 
oy Ox 
the equation of continuity is automatically satisfied. From 


/ 


y= V0p' xF(), = yey (8) 
Vy 
if follows that 
u = p'xF'(n), v= — V vp’ F(n). (9) 
Using (5) and (9), we have from (2) 
F" + FF" =F’ —14CF’, (10) 
where 
2 
==. au) 
pp 
The boundary conditions are 
u=v=0 at y=0; u= fx at y=oo. (12) 
These give from (9): 
F(0) = 0, F’(0) = 0; F'"(co) = 6/6’. (13) 
Using (6) and (11), we have 
2 
(By reba 
B B 
and the solution (positive) of this equation is given by 
Bas 2 
pe Cae ae (14) 
B 2 
_ Thus we have to solve (10) subject to 
| —C+VC244 
PAO} 0, ee 10) 0, (99) ; A A Fy) 


The equation (10) is non-linear and differs from that derived by 
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Blasius by the presence of the term CF’ representing the influence 
of the magnetic field. This equation subject to (15) has been inte- 
erated with the help of the Analog Computer EASE (at the Indian 
Institute of Technology, Kharagpur) for three different values of C, 


lOF 


i 


Fig. 1. Variations of F, F’ and F” with ». 


viz., 0, 2 and 3. The solution for the case C = 0 (i.e. the non- 
magnetic case) shows excellent agreement with Howarth’s so- 
lution. Graphs for F, F’ and F” are drawn in fig. 1 showing the 
influence of the magnetic field. From the graph it can be seen that 
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F’, and hence w/x, is essentially constant except in a layer of 
constant thickness. The constancy of this thickness is due to the 
decreasing pressure along the plate in the direction of flow, which 
compensates for its normal growth due to viscosity. It is also seen 
that f”’ decreases with the increase in the magnetic field. 

Using (9) the skin-friction 


( a 
Roy = | 
oy /y=0 


at a station x is given by 
V0 Ty 
xup’s!2 
The values of F”(0) for C = 0, 2 and 3 (as read from the graph) are 
shown in table I. 


eel) (16) 


TABLE I 
| F’’(0) 
1.25 


0.71 
0.58 


Cc 
) 
2 
3 


This shows that the effect of an increase in the magnetic field is to 
decrease the drag-coefficient. 
Again, from (9) the second equation of motion (3) becomes 


ov 1 op e2v 
VS i y ——_, 
oy poy dy? 
which can be integrated to give 
oe 1y2 a 17 
T= 90) de tos, (17) 


where ¢(x) is an arbitrary function. Using (5), the above equation 
gives 


$'(x) = — Bx. 
Thus (17) becomes 
jel uid) aes Pepe 
p Z 0 
which by virtue of (9) becomes 
7 12, 2 y 
ae aE i) ap Re). (18) 


p 2, 2 
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Since F(y) and F’(y) are known from the graph, the pressure 
distribution at a station « is known from (18). 

In conclusion I express my sincere thanks to Mr. A. K. Chaud- 
huri for his kind help in the solution of the differential equation 


with the analog computer. 


Received 5th July, 1960. 
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COMPLEX CONDUCTIVITY OF SOME PLASMAS 
AND SEMICONDUCTORS 


by Po H. FANG *) 


National Bureau of Standards, Washington, D.C., U.S.A. 


Summary 


The complex conductivities of plasmas and semiconductors have been 
calculated for several cases where the collision frequency can be expressed as 
a power function of the energy. From the result, some characteristic para- 
meters of the plasma originally investigated by Spitzer are estimated. The 
problem of determining the relaxation time from a non-symmetrical disper- 
sion is discussed. 


§ 1. Introduction. The nature of the electrical conductivity of 
weakly ionized gases and of some semiconductors is determined by 
the scattering processes. It is possible to express the effect of the 
scattering in terms of a relaxation time (or a collision frequency). 
Theoretically, the frequency spectrum of the a.c. conductivity can 
be derived from the energy dependence of the relaxation time 1)?). 
Hence, in principle, an experimental determination of the fre- 
quency spectrum should be capable of revealing the scattering 
mechanisms. Because of the presence of certain complicated 
integrals in the mathematical formulation, most of the theoretical 
discussions have been restricted to qualitative or asymptotic 
behaviour of the conductivity at extremely high or low frequencies. 
Consequently, the applicability of the theory in analyzing ex- 
perimental results has not been thoroughly investigated. In this 
paper the frequency spectra of the conductivity of several different 
kinds of gaseous plasmas, as well as conductivity due to various 
scattering mechanisms in semiconductors, will be calculated. From 
the results we will establish certain quantitative relations between 
the dispersion quantities and the energy dependence of the collision 


*) At present at the Institut Fourier, Grenoble, France, 
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frequency or the relaxtion time. Through interpolation we will also 
discuss the a.c. conductivities of the plasmas investigated by 
Cohen, Spitzer and Rotly, and by Spitzer and Harm °). 


§ 2. Mathematical formulation. The complex conductivity, o*, 
when the energy distribution is Maxwellian, can, in general, be ex- 
pressed as follows 1): 

3/2e—€ 
oo) = nee e3/ e ae (1) 

3o/am J v+io 
where is the number density of electrons, m and e are the mass 
and the charge of the electron, respectively, w is the frequency of 
the applied field, » is the collision frequency of electrons, and 
e = mv2/2kT where v, in the case of ionized plasma, is the relative 
velocity of electrons with respect to ions or molecules. The integra- 
tion limits here, and throughout this paper, are always from 0 to co 

for the corresponding integration variable. 

In the case of semiconductors, the relaxation time, which is the 
reciprocal of the collision frequency, is more frequently used be- 
cause of the usual quantum mechanical formalism. The expression 
is now 2) 


4e2% TE3/2e-€ 
—— 34/m m* 1 + tar 


o*() (2) 
where m in (1) becomes m”*, the effective mass of the electron. The 
relative velocity which is related to ¢ is now the relative velocity of 
the electrons with respect to the lattices or impurity centres. 


TABLE I 
Classification of Plasmas and Semiconductors according 
to p values 
pb Plasmas Semiconductors 
— 5/2 *) ref. 8) 
— 2*) ref. 8) 
— 3/2 ref. 3) LEIS) =25)4) 
—1 ref. 4), H20, NHs 
— 1/2 ref, 9) 
0) : Lorentz Lorentz 
1/2 ref. 5), He : ret) poo! 
1 ref. 6), air, ref. 18 7), No 
3/2 


*) The calculations were not made in these cases. See however table Ame 
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On the basis of experimental or theoretical results, various 
authors (see table I) have suggested the following relation between 
vy and e: 

Y= Ap leP (3) 
or, equivalently, 
T= Gye, (4) 


where the coefficient a» is independent of e, but may depend on the 
temperature. These relations are obtained either as asymptotic 
forms from the theory or, in the case of experiment, as approxi- 
mations over a finite energy range. In the integration, however, the 
Maxwellian factor e3/2e~* converges strongly at large values of «. 
At small values of ¢, the factor converges as the 3/2-th power of e. 
Therefore, the exact form of » is not important as long as » does 
not diverge more strongly than the (— 3/2)-th power of ¢ as ¢ ap- 
proaches zero. This point will be discussed further in § 6. Table I lists 
the values of # which correspond to plasmas of different gases and 
to semiconductors. In the case of semiconductors, specific scattering 
mechanisms have not been identified because several different 
mechanisms could lead to the same value of #. For the actual 
applications, therefore, the references should be consulted. The case 
corresponding to = O, which leads to Lorentzian dispersion, is 
well known, and the numerical evaluation is elementary. Two 
other cases, i.e., = +4, have been previously calculated 19). 
Also included in the table are cases of several values of which we 
cannot now relate to physical processes but these results might be 
useful in the analysis of the experimental results in study of plasmas 
and semiconductors. 


§ 3. Numerical evaluations. The numerical evaluation of the 
integrals of (1) or (2) is very simple when certain functions intro- 
duced by Dingle et al. 14) are used. These functions are: 


p(x) = (pt) | ae, (5a 
C(x) = (61) | 5 ae, (62) 


ePe-& 


Gy(x) = I | FF ae (Se 
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Substituting the relation of (3) into (1), 


CG = Kidp | a s <a eslees ede. (6) 
e2P + x2 

where K = 4e2n/3,/am, %p = apw. The integral of (6) may be 
represented by one of the functions of (5), (see Appendix 1). Nu- 
merical values may then be obtained from the tables of these 
functions by simple multiplication. In the extreme frequency 
ranges where the functions are not covered by the numerical 
tables, asymptotic forms, also given by Dingle, can be used. 
These forms may be simplified for special values of p (see Appendix 
1). The numerical values themselves are not listed since the calcu- 
lations involved are straightforward. 


§ 4. Analysis. The usual method of analyzing the frequency 
spectrum of a relaxation dispersion is to plot o’ and o”, the real 
and imaginary parts of o*, in the log w-space and to discuss the 
characteristics of the curves in this space. Fig. | shows the frequency 


r 4 
ol 1 es Sl 
‘ae tent aet 2 


Fig. 1. The frequency spectra of o’ and o”. 


spectra for various values of # in this representation. The half width, 
I’, and the maximum value of o’’ as a function of p are shown in 
fig. 2. The half width shows a minimum where p = O, and increases 
when the absolute value of # increases. Obviously, the two para- 


CONDUCTIVITY OF PLASMAS AND SEMICONDUCTORS Do) 


meters /’ and omax”’ alone are not sufficient to describe completely 
the line shape. For a comparatively more precise description, we 
introduced previously an Argand diagram 19), This diagram is the 
frequency trace in the complex plane of o*. We have shown that 
when # = 0, the trace of the diagram is a semi-circle. When 
p= +} the trace is approximately identical to the so-called 
Cole-Cole (abbreviated as C-C hereafter) arc 12). This arc is a 


i I pes 


2.0 (a) 


Fig. 2. The half width 7 and the maximum of o” for different values of p. 


semicircle with the center depressed from the o’-axis by an angle 
§ = x/2(1 — «), where « lies between O and 1. When the absolute 
value of p becomes larger than + 4, the arc approximation becomes 
poorer. Nevertheless if we draw arcs passing through the end points 
and the maximum point of the numerically calculated curves, these 
arcs provide a method for the description of dispersion data. When 
p is positive, the deviation from the C-C arc is greatest in the region 
of moderately high frequencies; it is greatest in the moderately 
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low-frequency region when # is negative. These diagrams are 
illustrated in fig. 3. Fig. 4 shows the dependence of 6 on #. 
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Fig. 3. Argand diagram of the complex conductivity: solid line, calculated 
from (6), dotted line. Cole and Cole arc approximation. 
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Fig. 4. The angle of Cole and Cole arc, 0, for different values of Bp. 
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The above description must be modified for extreme frequencies 
where o” approaches zero. In the case of the are given by (6), the 
tangents of the arc at w = 0 and w -> co make + 90° intersections 
with o’-axis (see Appendix 2). On the other hand, the tangents of 
the C-C arc make angles of + 2/2(1 — «). This difference will 
become important when considering the pressure dependence of 
the conductivity 1°), or the frequency dependence of the dissipation 
term in the Langevin equation 14), The terms in question diverge 
to infinity at infinite pressure or infinite frequency, respectively, 
for the C-C form of dispersion, while finite for the dispersion of (6). 


§5. Fully ionized plasma. Spitzer’s investigation of a fully 
ionized plasma has been much discussed. More recently, Sodha 
and Varshnil%) have given a simple approximation to replace 
Spitzer’s complicated relation for the relaxation time and the 
energy of the electrons. This approximation is for the case of the 
mean ionic charge Z = | and can be written as 
7 = Be0-722, (7) 
where B is a constant dependent on the impurity concentration. 
We have not explicitly calculated o* for this case because o* cannot 
simply be expressed in terms of known functions. However, it can 
be shown that, in general, o* is a continuous, non-singular function 
of p. Therefore, one can obtain the characteristic parameter of the 
complex conductivity by interpolation. The results are: 
o” max = 0.420 corresponding to r = 0.810¢9-722, 
0 = 8.4°, (8) 
I’ = 1.38 logarithmic cycle. 

The deviation of the arc in the Argand diagram from the C-C arc 


will be appreciable in the high frequency region. 
Spitzer has given the following relation for the case in which 


the ionic charge approaches infinity: 

Del, (9) 
The complex conductivity in this case has already been given in 
our discussion for p = — 1.5. 


§ 6. Relaxation time. In this section the relaxation time accord- 
ing to the usual formalism of relaxation theory will be discussed. A 
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dispersion quantity Q*(w) can be represented, in general,.by the 
following integral 14), 


2"(w) = 100) — OeeIn | HS — ar, (10) 


where g(r) is the distribution function of relaxation times, and N isa 
normalization constant such that 


N fe(r)dr = 1. (11) 


A similar formalism can be introduced in the investigation of the 
dispersion of the conductivity of plasmas or semiconductors 1°). 
Instead of dr in (9) or (10), sometimes d log 7 has been used in the 
integration in the literature of dielectric relaxation theory 1”). 
This difference has been considered purely as conventional and in 
our earlier investigation 1°), a formulation similar to that of (10) 
has been adapted. A unique formulation is provided by the di- 
mension of Q and g, and the requirement that the distribution 
function should be normalizable independent of the parameter # 
of (4) 15). In this way, a distribution function of relaxation times, 
g(r), is defined in the following manner 16): 


= o) = KN | EO —d log r. (12) 


Two applications of this relaxation time formulation are: 7) since 
the distribution function g(r) is the Stieltjes’ transform of o*(w), one 
can determine the distribution characteristics of relaxation times 
from the frequency dispersion data; and ii) from g(z), the average 
value of 7 or » can be evaluated. For the first application, in the 
present paper, we will only discuss the determination of a relaxation 
time parameter because the problem to find g(7) through the inverse 
transform is not required: g(r) is already given by (1), namely, 


sot (22)" of -()"], 


where (4) has been used, || stands for the absolute value of p. This 
additional notation, |p|, is used in order that the integration limits 
will always be taken from 0 to oo. A relaxation time parameter 79 
which corresponds to the most probable value of 7, is obtained ay 
maximizing g(r) with respect to log 7, 


70 = (5) ap. (14) 
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Now, experimentally, the relaxation time is usually determined 
from the reciprocal of the frequency corresponding to either the 
maximum of o’’(w), or the inflection point of o’(w). Denote the 
relaxation times determined in this way by 71 and 72, respectively. 
It can be shown !”) that the value of 79 is identical to that of 7 or 
72 tf and only if g(r) is a symmetrical dispersion function, in the 
following sense: 

g(z/t0) = g(TOo/T). (15) 


The distribution function of the present problem, which was given 
by (13), does not satisfy (15); therefore, 79 is not equal to 7, or 79. 
Table II lists the numerical values of to, 71, and re for different 
values of p. 

The distribution function g(7), as given by (13) has all the 
requirements of a distribution function, namely: (i), g(z) is a non- 
negative everywhere, (ii) it can be normalized, with a normalization 
factor determined from (13) and (15): 


N = 3/a, (16) 


and (iii), g(r) is zero both when 7 is at zero or infinity. All of the 
above characteristics of g(r) are independent of the value of sign 
of ~. In spite of this satisfactory behaviour, we observe that for 
certain values of # the mean value (first moment) of 7 or v does not 
exist. The mean value of 7, denoted by <7>g, is defined by 48), 


Ap fen eo des (17) 


leo 


<Tg =| rela) log 7 / fea log7 = ae 
The pathological nature of <7>g is obvious when one compares (17) 
with the Gamma-function integral and recalls that when 3— p < —1, 
the integral is not defined. One can apply an analytic continuation 
and define this integral in the sense of the Cauchy residue theorem, 
but the resultant <7>g has negative or infinite values for some 
p > 5/2. Therefore, the first moment of + remains unphysical. 

Molmud !4) has discussed a mean value of v, <y>m, which is 
defined as an average over the Maxwellian distribution; in the 
present notation, 


i 
uM = = |x em de, (18) 


a similar pathological problem as that of <7>g arises when p < — 3 : 
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The above discussion is dependent on the value of p. Now, p is 
dictated by the physical result, but otherwise arbitrary from the 
point of view of general mathematical formalism. Therefore, for 
some large values of the absolute values of p, the concept of the 
mean value of 7 or vy is not meaningful in the present problem 1%). 

As a numerical summary, a list of values of the various relaxation 
times as defined above and corresponding to different values of # is 
compiled in table IT. 


TABLE II 


Relaxation times calculated from different definitions 


p | T0/Ap 272 272 <T>g/ap 1/<v> Map 
— 1.5 S195 6.67 7.09 4.52 | 6) 
—1 2.50 3.08 3.00 2.50 | 0.50 
— 0.5 1.58 1.60 1.60 lol 0.89 
0) 1.00 1.00 1.00 1.00 1.00 
0.5 0.63 0.81 0.76 0.75 0.89 
1 0.40 0.88 0.83 0.67 0.67 
1.5 0.25 1.83 1.34 0.75 0.44 
2 0.16 = = 133 0.27 
25 0.10 — _— co 0.07 


§ 7. Conclusion. The complex conductivities of plasmas and 
semiconductors have been calculated in cases where the collision 
frequency can be expressed as a power function of the energy. 
From the results, we have established several numerical relations 
between the nature of the collision frequency and the characteristics 
of the dispersion functions. As an application of the results, we have 
estimated some characteristic parameters of the plasmas originally 
investigated by Spitzer. 

The problem of the determination of the relaxation time para- 
meter for the dispersion with non-symmetrical distribution of 
relaxation times, as well as the nature of the average relaxation 
time and the average collision frequency are discussed. 

The present mathematical procedure and the numerical results 
can be trivially extended to the calculation of the mobility in 
semiconductors due to the magnetic field, and to galvanomagnetic 
effects 2°) in general. 
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APPENDIX 1 


In the following, the complex conductivity will be expressed in 
terms of special functions of Dingle et al. The asymptotic forms 
of these functions will also be listed. o* will be expressed in the 
form: 


o* = KAap(o' — io”). 


The coefficient A, which depends on #, will be chosen such that 
o = | at the zero frequency. To avoid cumbersome subscripts in 
the discussion under special values of #, the subscript of x will be 
omitted. Also, the reciprocal of x will be denoted by y to further 
simplify the mathematical expressions. 


a) p= —3: 
A=3I, of =Gy(x2), 0” = x Case? 
When x is very small, 
3(x2) = 1 — 120% + 6.048 x 104x2, 
€4.5(x2) = 1 — 268% + 2.273 x 105x?. 
When «x is very large, we will expand €, in terms of y 
€3(x2) = y[0.167 — 0.201(y* — vy) + 0.004y — 0.018y log y], 
4.5(x2) = 0.025y(1 — 2.67y) — 0.0113y"/*exp(y!) 


[ — 2 exp(— 3+) cos (2 — * y)| 


b) p= —1: 
A =2.5!, Wass y2Co.5(y), o” = ty€3.5(y). 
For small values of y, 
Go.5(v) = x55 (1 + 4y?) — 0.945y2-8sin(0.785 ++ y), 
-Ca.s(y) = gig (1 — 0.335y2) — 0.220y?-5 cos(0.785 + 9). 
For large values of y, 
Co.5(y=) x2(1 — 15.75x?), 
€3.5(y) = «2(1 — 24.75x?). 
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4 =2, of =yPMhely 2), a = F yAe.5(y?). 
For small values of y, 
Wa.s(y2) = gy (1 — 0.667? + 1.333y4) — 0.9475 exp(y?), 
Wa(y2) = (1 — y2) — yt —.(0.577 + 2 log y)(1 — 94). 
For large values of y, 
Wo(y2) = x2(1 — 3x2 4+ 12x74), 
We.5(y2) = x2(1 — 3.542 + 15.754). 


A=15!, o =1/(1 + 2), 0” = x/(1 + 2). 
\p=+h: 
A=1, of = 2%o(x2), 0 = (1.5!)x 1.5 (x2). 


The asymptotic form of (x2) was given in c). For Y1.5(x?), 
when % is small 


W1.5 (4?) = s (1 — 2x2 — 4x4) + 1.965x8 exp(x?) ; 
when ~ is large 
Wi.5(e2) = y2(1 — 2.5y2 + 8.75 y4). 
fips ti: 
A =(—4)!, o! =F Cas(x), 0” = § x61.5(1). 


The asymptotic form for ©2,5 was given in b). For ©1,5, when x is 
small 


Ci.s(x) = a (1 — 0.667x2 + 0.1524) — 2.363x1-5 cos(0.785 + x); 
when ~ is large . 
C1.5(¥) = y2(1 — 8.75y2). 
2) pee 15s 
A =1, o' = (3!)y2G3(y2), o’” = (1.5)!yG1,5(y2). 
The asymptotic form of €3 was given in a). For small values of y, 
€1.5(y?2) = —2.667y2 (1 +0.761y2) +0.787y! exp(y')cos(1.047 —0.865y2). 
For large values of y, 
E1.5(y2) = 1: 39.4x2, 
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APPENDIX 2 


To show that o” makes + 90° intersections with the o’-axis, we 


observe, 
do” 7) | xete-? 5 lo eit De-€ 
=F : a ue de. 
om Ox J 62D +4 x2 | ex e2P +4 x2 


The numerator becomes 


0 gie-? 
(1 +x — ) 7 de. 
Ox e*P +. 2 


When x +0, the first term is finite, the second term and the 
denominator term approach zero from the negative value. There- 
fore, 


When x oo, all integrals become power functions of x with 
power — 2, and 
Co” 


do’ 


O(x) + 00. 


Therefore, o” intersects the o’-axis at + 90°, when w -> oo and at 
— 90°, when w = 0. 


Received 30th July, 1960. 
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ELECTROMAGNETIC GENERATION OF VORTICITY 
IN THE UNIFORM EFFLUX OF A CONDUCTING 
FLUID FROM THE SURFACE OF A MAGNETIZED 
SPHERE 


by J. D. MURRAY 


Mathematics Department, University College London, England 


Summary 

The electromagnetic forces in the flow of an electrically conducting fluid 
in the presence of a magnetic field are non-conservative and therefore 
produce vorticity. The simple case of uniform efflux of a conducting fluid 
from the surface of a magnetized sphere is studied. Two methods are de- 
veloped: one gives the solution for small conductivity and any magnetic 
intensity, the other gives the solution for small magnetic intensity and any 
conductivity. The case evaluated in detail is that for a magnetic dipole 
situated at the centre of the sphere. The magnetic lines, streamlines and 
vortex lines are found in closed form and shown in figs. 1 and 2 for two 
values of the parameters involved. 


§ 1. Introduction. In magnetohydrodynamics, the electromagnetic 
force is non-conservative and therefore produces vorticity. In all 
such flows of a conducting fluid the generation of this vorticity 
is of special interest. 

In 1), Ludford and Murray investigated the flow of an inviscid, 
incompressible, conducting fluid past a magnetized sphere: first 
order effects of the magnetic field and conductivity were studied. 
Due to the complexity of the solution for the vorticity, its behaviour 
was studied in detail only on the sphere, where it was found to be 
logarithmically infinite. The magnetic distribution was an axially 
symmetric one. In 2), Murray and Chi showed that in the two 
dimensional case the appearence of such a singularity in the 
vorticity was critically dependent on the magnetic configuration 
in the body. 

In 1) and 2), where the magnetic field originated in the body, it 
was shown that non-singular perturbations could not be made for 


ay Agee 
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small Ry, the magnetic Reynolds number. The non-uniformity 
at infinity could be incorporated in an exponential factor similar 
to that found in Oseen’s approximation in the viscous flow past 
a sphere. In the magnetohydrodynamic viscous flow past a semi- 
infinite flat plate Greenspan and Carrier%) and Murray *) 
show that an expansion in Ry is not possible when the magnetic 
field is externally applied in the direction of fluid flow at infinity. 

The investigation below is concerned with the steady efflux of 
an incompressible inviscid conducting fluid from the surface of a 
sphere in which a magnetic field originates. The methods are de- 
veloped for a general magnetic distribution. In this situation, for 
zero conductivity, the fluid velocity and magnetic intensity tend 
to zero algebraically at large distances from the sphere. The 
simplest case in which the electromagnetic body force is non-zero 
is that in which the fluid field (for zero Ry) is represented mathe- 
matically by a source at the centre of the sphere, which is con- 
sidered to be uniformly magnetized, i.e. a magnetic dipole is 
situated at the centre. Here the undisturbed fluid velocity vector 
is not parallel to the magnetic intensity vector and vorticity will 
thus be created. For such a configuration the current lines and 
vortex lines are circles. 

In the case discussed in §3 below, it is shown that a simple 
perturbation in the magnetic Reynolds number Ry is uniformly 
valid for all values of 8, the ratio of a representative magnetic 
pressure to a representative fluid dynamic pressure. On the other 
hand, for small values of 8, a simple perturbation can be found 
in B for all Ry (see § 4). If the latter solution to first order in f is 
then expanded for small Ry, it agrees with the former. Closed form 
solutions can be found to any order of Ry or f for the magnetic 
field, the electric current, the fluid velocity and the vorticity. The 
case worked out in detail is that for small Ray: no difficulties are 
encountered for large Ry. The streamlines, constant magnetic field 
lines and vortex lines are shown in figs. 1 and 2 for various values 
of the parameters involved. 

The flow is completely independent of the conductivity of the 
sphere as long asit is finite. If the magnetic field is originally assumed 
to be frozen into the sphere (i.e. the sphere has infinite conductivity), 
the results are different: these are given in §'5. 

Astronomically it is of interest to consider pressure balanced 
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magnetic fields in the atmosphere of the sun (see e.g. ®)), which 
may be considered to be a uniformly magnetized sphere. In this, 
the efflux of an infinitely conducting gas from the surface distorts 
the magnetic lines and pushes them outwards from the centre like 
elastic strings under tension, until the magnetic pressure balances 
the fluid pressure. If the conductivity of the gas is finite, the field 
lines are still pushed out, but the steady state then consists of the 
fluid passing across the field lines in their displaced position. This 
is similar to the problem discussed below. Accordingly, the case 
worked out in detail is that for a dipole situated at the centre of 
the sphere. With the method described below, the fluid pressure 
distribution on the surface of the sphere can be varied easily by 
addition of further fluid sources and sinks inside the sphere as long 
as axial symmetry is preserved. 


§ 2. Basic equations and their non-dimensional form. Neglecting 
the displacement current, the equations governing the steady 
flow of an electrically conducting, incompressible, inviscid fluid 
of constant properties are 


div q = 0, (1) 

curl g x q = — grad (p/p + 492) + a < H, (2) 
j = curl H = o(E + xq X 4), (3) 

div H —0 = ‘curl E, (4) 


where yw and o are the magnetic permeability and conductivity 
of the fluid respectively, and q, p, p, jf, H and E have the usual 
meanings. 

The specific three-dimensional geometry considered is that of a 
sphere of radius a, situated at the origin y = 0, where 7 and @ are 
the axially symmetric polar coordinates as in fig. 1. The uniform 
efflux of fluid is simulated mathematically by a source at the 
_ origin. The magnetic field H, which originates in the sphere, and 
the velocity q lie in the meridian plane and are independent of the 
azimuthal angle (see figs. 1 and 2). From the conduction equation 
(3) the electric field E is perpendicular to the meridian plane, and 
from the second of (4), E = 0 if it is to be finite on the axis. 

The quantities g, H and 7 are made dimensionless by referring 
them to the velocity U at the surface of the sphere, a representative 
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magnetic intensity #, and the radius of the sphere a, respectively. 
The equations in non-dimensional form become (retaining the same 
symbols for their non-dimensional form) 


div q =0, (5) 

curlg X q = — grad P+ f curl H x A, (6) 
curl H= Ryq X A, (7) 

div H = 0, (8) 


where P = p + 4q?, B = pwh?/pU2, Ry = Uauo and ? is now the 
pressure divided by pU?. 


Fig. 1. Magnetic and fluid field lines when Rm = 0.2, B=5, w’ = 2y 
and o’ is finite. 


When Ry = 0, or 8B = O, the streamlines are straight lines and 
the velocity is given by 


d = Jo = (I/r*, 0, 0). (9) 
In the axially symmetric case, H = (H;, Hg, 0) and 
eee 1 @A 
Ay = — eet Hg aes se 
"sing a ° - ° yrsinO or ’ un) 


where A is a function of 7 and 6 alone. 
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. § 3. The perturbation in Ry. For small conductivity we expand 
in powers of Ry and write 
edo Mi. = Ay + Ry, = ...: 
ely . b=fot+ Rufic+... (11) 

Substitution of (11) in (5) to (8) gives go as in (9) and the non-zero 
component of (7) as 

027A sing 0 ( ] se 

cr? r 06 \sin6 & Ee 


(12) 


y=0 


/W=—342 


Fig. 2. Magnetic and fluid field lines when Rm = 0.2, B=5, pw’ = 2u 
and o’ is infinite. 

Inside the sphere E = 0 again, and so curl H = o’E = 0, where o’ 

is the conductivity of the sphere. Thus, inside the sphere Ao (and 

also A;, 7 > 1) satisfies (12). 

Any axially symmetric distribution which is a solution of (12) 
can be taken. In the specific case here considered, i.e. that of a 
dipole situated at the origin pointing in the direction 6 = 0 (figs. 1 
and 2), the potential Ao which satisfies the boundary conditions 
on the sphere (continuity of the normal component of magnetic 
induction and tangential component of magnetic field) and tends 
to zero at infinity is 


By Sine G4 LOT 72), 
Ao(r, 0) = 


a (13) 


sin? 6 — wu 3 
: M for sin? 0, fory < 1, 


Y PATE wT 


Ag, o).= 
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where yw’ is the magnetic permeability of the sphere. 
The equations for q; and Hy are 


div qi = div Hi = 0, (14) 
curl qi X qo = — grad P; + fcurl Hi x Ho, (15) 
curl Hy = qo X Ho, (16) 


H, and qi are found simultaneously. The non-zero component of 
(16) gives the equation for the potential A; (of M1) as 


027A; | sin o ( 1 a ( Su = (17) 
ov? ov?) 8 \sind §26-, oi 


a particular integral of which is 


2 ( Vis. ) sin? 9 
4\2u+ yu ye 


To this solution must be added the appropriate complementary 
functions of (17), (i.e. sin? 0P’,(cos 6)/7”, n an integer and P,,(cos 6) 
the Legendre polynomial), so that the boundary conditions on the 
sphere are satisfied to O(Ry). Inside the sphere, irrotational 
disturbances must be added to Ap of O(Ry) of the form 


Ryr” sin?6P’n(cos 8), 


which are solutions of (12). A little manipulation shows that the 
appropriate solutions satisfying the boundary conditions are 


3) ; sin? 6 [ 2(u x l 
Ay= ( La 7 [ areal Obs i: forr >1 
4\ 2u + u y (2u + uw’) Y 


1= % <3 lle —y#2sin29, forr <1 
2 (2u + w’)? Boye 


With Ao and A, from (13) and (18), A = Ao + RA; satisfies the 
equations and boundary conditions to O0(Ry) and from (10) the 
components of magnetic field are given by 


He=(5 Su’ jee (4 (u + w’) ==] 
2p + re (2u + wu’) Y 


Hy = (5 =) ti i al (uth) “| 
oe es eee PMI 2uek uh oo 


(18) 


(19) 


? 
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outside the sphere, and 


H 2cos6 | ( b— we ) eer Lie ye P 
= = COS / = —= SU 
i 73 2u + mw ers (2u + Fes : , 
sin 6 wu mw bee’ ad 
,= - —})2sin 6 — Ry 3 ———_—— sin 9, 
78 2u+ bu (2u + pw’)? 


inside the sphere: (18), (19), and (20) are correct to O(Ry), 
The vorticity vector, curl qi, is perpendicular to the meridian 
plane and has therefore only one component w . From (15) 


OP, 
v= — 7 ais BF, 
or 7 
OM) uy 1 OP, BF ( 
jae divin asic, oc. 
where F = (F;, Fg, 0) = curl Hi X Ho = (Go X Ho) X Ho, so 
Su’ 2 sin? 6 Su’ 2 sin 20 
a ( lu | sin Fy =( HM = ae 
2u+u hi Qu + bu v 


From (21), since , is finite at infinity, 
6 Si6 2 sin 26 l Su’ 2 sin? 0 
o = 62(5~—) Pr =8(5_) == @ 
7\2u +4 76 DN 2a ih if 
Thus, @ = (0,0, w) = (0,0, Ry) to O(Rm), with w, given by 
(22) for all 8. Lines of constant vorticity are shown in figs. 1 and 2. 
The velocity vector is easily found from qo and 1, where 1 1s the 
Stokes stream function for qi and 
Cy Sig. «a ( l Ow1 
or? 72 06 \sind 26 
with w, from the first of (22). The boundary conditions on yj are 


wy, > 0Oas7 — oo, and the normal component on the sphere is zero. 
The appropriate solution of (23) is 


Bu" 2 sin? @ cos 1] 
7\ 26 + pb 72 Y 


) = —yrsin $0, (23) 


and 


Qi = (M1y Y1o» 9) | 
a Hf Bu Me asap a i) sin Ze- a) o|. (25) 
rae) aay v4 v ors r 
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The streamlines are given by y = yo + Ruyi = constant, where 
yo = — cos 0, giving qo = (1/72, 0, 0) and are shown in figs. 1 and 
2 for two values of o’, and given Ry and f. Note that (24) and (25) 
hold for all p. 

Higher order terms in Ry can be easily found, also, for all . 

Note that the ratio of any term to the preceding one in the q 
series tends to zero as 7 — co, whereas the similar ratio in the H 
series tends to a term of O(Ry). This is not, however, surprising 
since in the interaction of Hp and qo, Hp = O (1/r3), and qo = O(1/r?) 
for large 7. 


§ 4. The perturbation in B. In the above section, closed form 
expressions were obtained for H, gq and @ for all 6 and small Ry. 
If, on the other hand, solutions are sought for all Ry, a simple 
expansion for 6 can be taken, similar to that in 1) and 2). Here we 
write 


q=d+ q+... H= Ho + pM + ..., dp = fo + Boi + ... (26) 


Here, Ho is not the undisturbed field since, to get the dimensional 
magnetic field, Hp must be multiplied by h, which is proportional 
to 6, and it is this quantity which would tend to zero as B +0. 

On substitut ng (26) in (5) and (6), qo is again given by (9), and 
(7) and (8) become 


curl Hp = Ruqo X Ho, div Ho = 0. (27) 
Using (10), the non-zero component of the first of (27) is 


OPAg _ sind 0 ( l a R 1 @Ao 
2 "60 sin (00. JY Fan Oe or 
simple product solutions of which are 
Ao(r, 0; Ru) = Fn(r/Rm) sin? 6 P’n(cos 8), (29) 


where m is a positive integer, P’» the derivative of the Legendre 
polynomial, and 


Fn(r/Ru) = (Ru/r)" fn(Ru/r) exp (— Ru/r), (30) 
with 
adit ol (1+ 2) Ru | (n+ 2)(n+3) /Ru\2 
LS la oe 1!(2n + 2) eran! mee 


FLOW VORTICITY IN THE FIELD OF A MAGNETIZED SPHERE 723 


(see Whittaker and Watson §), Chapter XVI, ). For fixed 7, 
Ry"Ao > sin*6P’y (cos @)/7™ as Ry +0; 


i.e. solutions of curl Ho = 0, as it must. 

The solution of (28) is now found by adding to the appropriate 
solution (29) complementary functions of (28) so that the boundary 
conditions on the sphere are satisfied for all Ry. For a dipole 
situated at the origin, the appropriate solution to take is 


sin2 6 


Av=K 


; fi(Raz/r) exp(—Rwm/7), (7 = Os (31) 


the solution x = 1, where K must be determined from the boundary 
conditions on the sphere. From (30), 


SREY, 3 ea 3 (ers 
sil eT | 306 
fi(Ru/7) is 4 oY | leks Y (PE Bh3 Vy | 


The potential inside the sphere is of the form 


ait 
Bet Pai! ye 1), (32) 
Y 


where & is determined in conjunction with K. The representative 
field h can now be identified as M/a3, where M is the moment of 
the dipole. Using (31) and (32), the continuity of the normal 
component of magnetic induction and the tangential component 
of magnetic field give 


3u’ exp (Ry) 
~ 2ufi(Ru) + w'{f(Ru) — Rulf(Rm) — fi’(Rm))} ’ 
(uw — w’) fh(Ru) + Ru’ [A(Rm) — h'(km)] 
~ Qui(Ru) + w'{fi(Ru) — Ralf(Ra) — fr'(Ra))} 


(33) 


where the dash on f; denotes differentiation with respect to the 
argument. The components of the magnetic field are given by 


Ho, = K : = 5 fi(Ru/r) exp (— Rm/r), 
Y 
| | (34) 
pi ,,— wnt (Ran) “Mef(Ru ft (Rui) | (exp(—Bai7), 
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outside the sphere, and by 


Y 6) sin 9 
Ho, = me NMA, daly. 
73 


— 2k sin 6, (35) 


inside the sphere, with K and k from (33). 
When Ry is small, 


fi(Rmu/7) exp(—Ry/r) = 1 — Ru/r + 0(Rmu2Z/r?), 


aie 
(cH) f+ x tt] + out 
Zea yh 2(2u + pw’) 


hale / 3R 1 nf 
=( lu Ve f al ] + O(Rm?), 
2h + | 4(2u + w')(u — b’) 


and (31) and (32) give the same result as the sum of (13) and (18) 
to the same order in Ry. The vorticity m1, the stream function 1 
and the velocity gi are then given by (22) (24) (25) respectively 
to 0(Rm), and the magnetic, fluid and vortex lines are as in figs. 
1 and 2, even though f is not small. 

For all Ry, the vorticity @1; = (0,0, w1), from (26) and (6), is 
the solution of 


@O1 X qo = — grad Py + Ru(qo x HA) x Ho, 


where Ho is now obtained from (31). The non-zero components 
of the last equation give 


OP, Ru ‘ asl oP, R 


Sl po ee oe ee ca! 
or y2 ( 00) , y 06 y (Hog H,), 


with Ho,, Hog from (34). The non-zero component of the vorticity 
@ is, to 0(A), 


@ = fo, = Rees si ae oa fi(Rarlr) - 


*(fi(Ra/r) — (Ru /”){f(Ru/r) — fr'(Ru/r)¥) + 


exp(—2Ry/r 
+7 ear) alta Rare) Ranier] (36) 


co 
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where the lower limit, infinity, in the integral ensures finiteness of 
@ as 7 — oo. For small Ry, (36) gives (22) to 0(Rm). For large Ruy 
the /;(m/r) series in (30) is used: the series converges quite rapidly. 
Equation (36) can be integrated term by term with /1(Rym/r) taken 
up to the last term used in the summation to obtain Ay. The stream 
function and velocity can be found in the same way as in § 3 to any 
degree of accuracy required for any finite Ry. 


§5. The perfectly conducting sphere. The above analysis is for 
a sphere of arbitrary conductivity o’, however large, as long as it 
remains finite. The final expressions for the magnetic and fluid field 
do not involve o’. For the case of infinite conductivity, the magnetic 
field is frozen into the sphere and the boundary conditions on the 
magnetic field are different as are the results. 

For the dipole case above, the field inside the sphere is always 
given by H = curl(0, 0, 7-1 sin29). The only condition on the field 
outside the sphere is now that of continuity of the normal component 
of magnetic induction: there is no condition on the tangential 
component of magnetic field. The first equation of (18), the first 
of (22) and (24) become respectively 


yee \esin® 6 l 
Ay=-—|— 1 ; 
4\ u fs ry 


6 / mw’ \? sin 26 
a: 52 lt ) yo” i 


2 [{ uw’ \2 sin? 6 cos 0 1 
CAGES) Y vy 


and (31) becomes 


ae EUs emir ee) 


r film) 
while the vorticity is the same as (36) with K? replaced by 


( a i exp(2Ran) 
w/ [fi(Ru)}? 


Thus, the assumption of infinite conductivity in the sphere would 
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lead to an error, in the magnetic field say, in § 4 by a factor 
K — 
(w'/u) fi(Rm) exp Kus 


Su 
~ f(Ru)[2ef(Ru) + #'{f(Ru) — Rul f(Ru) — fr'(Ru)]}]” 
and comparable factors in the fluid field. 
Figs. 1 and 2 illustrate the difference in the magnetic field and 
fluid field for Ry = 0.2, B = 5, w’ = 2, and o’ finite and infinite 
respectively. 
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Summary 


It is usual in the analysis of one-dimensional channel flows to study the 
behaviour of the analogous isentropic flow since, first, it retains the essential 
features of flows of practical interest and, secondly, it is simpler to describe. 
Although in conventional channel flows it is sufficient to neglect heat 
addition and friction to ensure isentropicity, in the MHD case it is in addition 
necessary to neglect Joule heating. This is accomplished by considering 
the fluid as having infinite electrical conductivity. However, this procedure 
does not necessarily imply infinite currents, since the external resistance will 
limit current flow. In the conventional problem, if we assume an isentropic 
flow, we are able to obtain a once integrated form of the governing equations. 
Such once integrated solutions are not possible in the present isentropic MHD 
channel flow, but equally simple solutions can be found and are presented. 
Examples of application of these results to the crossed field MHD generator 
and accelerator are also given. 


§ 1. Introduction. Application of the crossed field generator 
to the large scale generation ot electrical power 1) and of the crossed 
field accelerator to space propulsion 2) has stimulated considerable 
interest in MHD channel flows in general. The purpose of this paper 
is to propose a simple procedure for studying their behaviour. 

As is well known, a useful approximation for the analysis of 
channel flows is to assume that the flow is quasi-one dimensional. 
Such an approach permits a considerable simplification and will 
be followed in the present study. We also know that the analysis 
of the resulting one-dimensional equations will be complicated by 
the presence of the Lorentz force and Joule heating terms °). 
Therefore, a further simplification would be to neglect Joule heating. 
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If we also neglect heat conduction and viscous dissipation, and if 
we require that changes of state occur very slowly, the flow may 
be considered isentropic. . 

However, the Joule heating can only be neglected if the fluid 
conductivity is assumed to be infinite. Normally, in the study of 
unconfined flows this assumption implies no relative motion 
between the fluid and the magnetic lines of force, since any such 
motion would induce infinite currents. In the present confined 
channel flow, though, such relative motion is possible, since the 
current flow is limited by the external resistance. Therefore, in 
the following analysis we will make such an assumption and study 
the resulting, simplified equations. 


§ 2. Analysis. The basic equations for the steady flow of a. 
perfect, non-viscous, non-heat conducting gas of infinite electrical 
conductivity can be written as follows: 


momentum: 
pv Viv+tVp=t x B, (1) 
isentropicity : 
V(p/p’) = 9, (2) 
continuity: 
V-(pv) = 0, (3) 
Ohm’s law: 
E+vx< B=0, (4) 
Maxwell’s equation: 
i= —V x B, (S) 
Me 


where we have used the rationalized MKS system of units and 
where v is the gas velocity vector, p the gas pressure, p the mass 
density, y the isentropic constant, i the current density vector, E 
the electric field intensity vector, B the magnetic induction vector 
and we the magnetic permeability, which is assumed constant. If 
we recall that the total magnetic field can be considered as the sum 
of an applied and an induced field and observe that within the channel 
the applied field must be irrotational, then (5) can be written as 


l 
i=—V xd, 6 
Le (6) 


where b is the induced magnetic induction vector. 
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Next we consider the flow through a channel of slowly varying 
cross-sectional area, such as that shown in fig. 1. For this geometry 
we employ the usual assumptions of quasi-one-dimensional flow. 
That is, we consider that there is only an x-component of velocity 
u, a z-component of applied magnetic field B, y-components of 
electric field E and current density 7, and that these along with the 
pressure and density are functions only of *. The z-component of 
induced magnetic field has been shown by McCune and Sears 4) 
to be negligible for infinitely long channels of slowly varying height 
and, therefore, will be neglected. They have also pointed out that 
the x-component of induced magnetic field within the channel 
varies linearly with and is an odd function of z, so that its average 
value is zero. Since our one-dimensional approach considers only 
average quantities, this is consistent with our assumption of a 
pressure ~ which is a function of x alone. 


Magnetic Field 


Electrode 


Fig. 1. Configuration of general MHD channel with crossed electric and 
magnetic fields. 


If we then neglect both components of the induced magnetic field, 
we obtain the following simplified set of equations. 


pus + oe = 3B, (7) 
= (lpr) = 0, 6) 
& (puA) = 0, 9) 
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where A is the channel cross-sectional area, which can be a function 
of x, and (9) is the continuity equation appropriate for the study of 
one-dimensional channel flow. 

Before entering into a discussion of particular solutions to these 
equations, some general observations would be of value. The variable 
E can be omitted from consideration initially, since with B a given 
function of x, it can be found once w is obtained from the solution 
of (7), (8) and (9). There remain, therefore, three equations to be 
solved for the five variables uw, p, p, A and 7. Thus, in order to solve 
the problem, we must specify the x-dependence of two of these 
variables and give the initial values of the remaining three. 

Equations (7), (8) and (9) can now be rewritten as follows: 


du dp : 
pO ee (dal) 
POR 25 (12) 
po” py 
pouoAg = pua. (13) 


If the Lorentz force term 7B in (11) were to be taken to be zero. 
then (11), (12) and (13) would describe a conventional isentropic 
flow. In this case, the equations can be integrated once and become 


wi ym es " 
2 ey eee Pt Te : 
Po 
po” pY 
powuoAg = pu, (16) 


where the subscript zero denotes the value of the variable at some 
initial station. These equations can be handled in two different ways. 
First, we can specify the value of one of the variables at some down- 
stream station and solve for the value of the other three variables at 
that station from the above equations. In this procedure, the distance 
between the initial and downstream stations is not specified, 
although it cannot be zero if all variables are to vary slowly with x. 
Another procedure would be to specify the x-variation of one of the 
four variables and solve for the other three as functions of x from 
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(14), (15) and (16). In this case, the distance between the initial 
station and one downstream is clearly defined, and although the 
calculation is not quite as easy, it is considerably more informative. 

When the Lorentz force term is included to permit the study of 
isentropic one-dimensional MHD flows, a once integrated form of 
the equations, such as (14), (15) and (16), cannot be found. Solutions 
for the variables as functions of x can, however, be obtained by 
giving their initial values and specifying the x-variation of two of 
them. In the next section, two simple calculations will be made to 
illustrate this procedure. 


§3. Examples. As a first example, let us consider a linearly 
increasing velocity and area so that we know 
u—=—u+ax; A=Ao+ Ox. (17a, b) 


In this case, (13) yields the variation of density with respect to x: 


A 
pottoAo (18) 


ae abx? a (aAg + buo) x + ugAo ; 


Combining (12) and (18) gives an expression for the pressure which is 


przls polnoe ey . (19) 
[abx? + (aAo + buo) x + upAo]” 


To solve for 7, we introduce u(x) from (17a) along with p(x) and 
p(x) (equations (18) and (19)) into (11) and find the following result: 


~~ els} poloAo(uo + ax) a 
See? | abx2 + (aAo + buo) x + updo 
po(uoAo)” y(2abx + aAo + buo) (20) 


Ps [abx2 + (aAo + buo) x + uoAol”*t 


where B is a function of x which remains to be specified. The channel 
geometry as well as the power consumed or generated can be 
determined with the aid of (10). We first define the electrode height 
as h and the spacing between electrodes as g, so that we can write 


A, = Ph. (21) 
Then returning to (10), we substitute for u(x) from (17a) to obtain 


E'= Vig = (uo + at) B, (22) 
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where V is the constant voltage between the electrodes. The 
electrode gap is, therefore, given by 


V 


= (23) 
(uo + ax) B 
The electrode height is thus available from (21) and (23): 
rice (Ao + bx)(uo + ax) B (24) 


V 


Finally, the total power generated or consumed is equal to the 
voltage V times the current flowing. The net current flowing for 
a channel of length L is simply the integral over x of the product of 
current density 7 and electrode height h. Therefore, the net power is 
calculated to be 
a2 
1P == pouo2A oaL + pouoAo nae Sarat 
yPo(uoAo)” 
y — | 

This relation shows us that the net power is independent of the 
voltage V and applied magnetic field 6b and depends only on the 
initial conditions, the total channel length, and the assumed forms 
of wu and A. The field 6 cannot, however, be made arbitrarily small 
since, according to (20), this would cause the current density to 
become extremely large. Such extreme current densities cannot 
be obtained in practice due to the fact that electron emission from 
electrode surfaces is limited. The constant voltage V primarily 
influences the values of g and / and, if not chosen properly, can make 
their ratio impractically large or small. In particular, it cannot be 
zero if B ~0, since then g would be zero and h infinite. 

For our second example, we consider a case for which the current 
density 7 is given as a function of x. We then specify that the 
velocity wv is constant, the applied magnetic field B is constant and 


and that the current density 7 increases linearly with x«..We can 
then write 


[(uwoAo)1-” — (abL?2 + (aAo + duo) L + upAo)!-7). (25) 


1= 19 + cx, 


and since wu is constant, (11) reduces to 
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Integrating to obtain the pressure p we get 
b= po + Bx(to + cx). (27) 
The density is obtained from (12) and is given by 


0 


ee 
pe “poll? [Po + Bx(to + 3cx) |", (28) 
and the area A is found by combining (28) and (13): 
Aopol/Y 
AS P 4 (29) 


[po + Bx(io + $ex)}/7 
The power P can be calculated if we first obtain g from (10). 
Recalling that E = V/g, we have 
g = V[uob, (30) 


and the gap is constant for this case. Then we find / from (21), 
(29) and (30): 
oe Aopol/’upb 
[po + Bx(to + dcx)” V © 


(31) 


Finally, the net power is calculated to be 
A oPol/Yuoy 
ea | 


P= 


{| po 1 BL (to + $cL)|v-UyY) = poY—Uy)}. (32) 


We find in this case, in contrast to our previous result, that the 
net power P depends on the applied magnetic field B. This situation 
occurs because the power per unit duct volume is £7 and E varies 
directly as B (see (10)) while 7 is assumed to be independent of B. 
In the previous case, 7 was inversely proportional to B (see (20)), 
so that the product Ei was independent of B. 

Although the solution of the second case was quite simple, it 
should be pointed out that, in general, when 7 is one of the two 
variables specified, the equations will be difficult to solve. On the 
other hand, the first case was treated readily, and any other case 
in which we specify two of the four variables Pp, p, and A will 
also be solved easily. Since these latter cases are those which would 
be of most interest for practical applications, the concept of 
isentropic one-dimensional MHD channel flows developed here 
should prove useful. 


84 ISENTROPIC MAGNETOHYDROMAGNETIC CHANNEL FLOW 


§ 4. Conclusions. In the preceding sections, we have formulated 
the one-dimensional approximation to isentropic MHD channel 
flow and have shown how solutions to the resulting equations may 
be obtained. These solutions can be interpreted as describing either 
a crossed-field accelerator or generator depending on whether the 
net power P is negative or positive. For example, if in the first case 
considered in the last section, where velocity increases linearly with 
x, we select our constants such that P is negative (power added to 
gas), then this is an example of a useful crossed-field accelerator. 
If P were positive, the device would be a generator of rather 
unlikely design. The second case could, on the other hand, be con- 
sidered a generator (since velocity is constant) by taking 79 to be 
some negative value and limiting the duct length L to a value less 
than the length at which the current density 7 is reduced to zero. 

In either application, the present approximation will be valuable 
for the description of high efficiency devices. 


Received 19th September, 1960. 
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